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ABSTRACT 

In  order  to  study  the  nucleate  pool-boiling  performance 
of  R-114  and  R-114-oil  mixtures  from  enhanced  evaporator 
tube  surfaces,  an  experimental  apparatus  was  designed, 
constructed  and  instrumented.  The  evaporator  was  made  of  a 
T-shaped  Pyrex  glass  container.  Boiling  occurred  from  a 
smooth,  hard-copper  tube,  15.9  mm  (5/8  in.)  in  outer  diam- 
eter, 12.7  mm  (1/2  in.)  in  inside  diameter  and  431.8  mm  (17 
in.)  in  length.  The  tube  was  heated  using  a  cartridge 
heater,  and  it  was  instrumented  with  8  thermocouples  to 
measure  the  wall  temperature.  A  Hewlett-Packard  3497A  data 
acquisition/ control  unit  and  a  9826A  computer  were  used  to 
collect  and  process  data.  The  condenser  was  cooled  by  an 
ethylene  glycol-water  mixture,  which  was  maintained  at  about 
-17  °C  by  means  of  an  R- 12  refrigeration  system.  Nine  data 
runs  were  completed  to  de-bug  the  experimental  apparatus  and 
to  check  for  reproducibility.  During  .all  data  runs,  espe- 
cially at  higher  heat  fluxes  (greater  than  10  kW/m2 ) ,  large 
temperature  variations  were  observed  along  and  around  the 
active  boiling  length  of  the  test  tube.  The  data  were 
compared  with  data  found  in  the  literature  and  reasonable 
agreement  was  obtained. 
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NOMENCLATURE 

A  area 

A  tube  outside  surface  area  of  active  boiling  section 

A^  cross- sectional  area  of  tube 

A  tube  outside  surface  area  of  non-boilng  section 

ag.  .  .  Brj  conversion  coefficients  for  thermocouple 

CL  specific  heat 

D  diameter 

D-  tube  inside  diameter 

1^  tube  outside  diameter 

Dl  diameter  at  the  position  of  the  thermocouple 

D>  outer  diameter  of  the  boiling  tube 

DCP  discrepancy 

E  thermocouple  reading 

g  gravitational  acceleration 

h  heat- transfer  coefficient 

hf  latent  heat  of  vaporization 

I  current 

I  current  reading  by  AC  Current  Sensor 

k  thermal  conductivity  of  liquid 

k  thermal  conductivity  of  copper 

L  active  boiling  tube  length 

L  non-boiling  length  of  the  test  tube 

Nu  Nusselt  number 

P  pressure 

P  critical  pressure 
c  r 

p  tube  outside  wall  perimeter 

Pr  Prandtl  number 

Q  heat- transfer  rate  from  boiling  surface 

Q  heat- transfer  rate  through  one  unboiling  end 

Q  heat-transfer  rate  from  cartridge  heater 

H 

Ra  Rayleigh  number 
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T  temperature 

T  average  wall  temperature  at  the  thermocouple 

location 

T\  calibration  bath  temperature 

T  temperature  of  the  thermocouple  location 

Tsat  saturation  temperature  of  the  boiling  liquid 

T"  outer  wall  temperature  of  the  boiling  test  tube 

V  voltage  across  the  cartridge  heater 

V  voltage  reading  by  AC-DC  true  RMS  converter 
a  thermal  diffusivity 

3  volumetric  thermal  expansion  coefficient 

6  uncertainty  in  measurement  and  calibration 

y  viscosity 

v  kinematic  viscosity 

$  volume  fraction  of  pure  refrigerant 

pL  saturated  liquid  density 

p  saturated  vapor  density 
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I.  INTRODUCTION 

A.  BACKGROUND 

In  order  to  reduce  the  size  and  weight  of  air- 
conditioning  systems  and  water  chillers,  the  United  States 
Navy  recently  has  considered  two  different  approaches:  1)  to 
change  the  type  of  the  refrigerant  from  R-ll  to  R-114,  and 
2)  to  use  enhanced  evaporator  tube  surfaces,  such  as  Gewa-T, 
Thermoexcel-E ,  and  High-Flux.  Schematic  drawings  and  photo- 
micrographs of  these  enhanced  surfaces  at  various  magnifica- 
tions are  given  in  Figures  1.1  ,  1.2  and  1.3.  A  combination 
of  R-114  and  enhanced  surfaces  may  lead  to  significant 
improvements  in  system  efficiency  as  well  as  to  smaller  and 
lighter  components  in  future  Navy  applications. 

B.  ADVANTAGES  OF  USING  R-114 

The  latent  heat  per  unit  volume  of  R-114  is  50  BTU/ft3 
at  70  °F.  As  a  result,  the  heat  carried  by  each  cubic  foot 
of  R-114  is  nearly  twice  as  much  as  that  carried  by  R-ll 
[Ref.  1]  at  the  same  temperature.  In  other  words,  the 
volume  of  R-114  vapor  circulated  per  ton  of  refrigeration  is 
about  half  that  for  R-ll.  Inherently,  the  size  of  the  system 
components  can  be  reduced  significantly  using  R-114  rather 
than  R-ll.  R-12  and  R-22  in  fact  can  transfer  the  heat  more 
efficiently  than  R-114,  but  they  require  higher  pressure 
ranges.  However,  operating  pressures  for  R-114  are  moderate, 
similar  to  R-ll,  such  that  light-weight  components  may  be 
desirable  and  less-experienced  personnel  may  be  involved  in 
the  installation  and  maintenance.  While  the  refrigerants  in 
general  are  nonflammable,  R-114  and  R-12  are  rated  to  have 
the  lowest  toxicity  among  all  refrigerants  [Ref.  2].    Based 
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on  the  above  discussion,  R-114  is  probably  the  most  suitable 
refrigerant  for  use  on  board  submarines  and  ships . 

Further,  R-114  is  considerably  more  stable  toward  heat 
than  R- 12 .  All  standard  ferrous  and  nonferrous  metals, 
except  zinc  and  magnesium,  and  aluminum  alloys  containing 
appreciable  amounts  of  zinc  and  magnesium,  may  be  used  with 
R-114  [Ref.  3].  The  physical  properties  of  R-114  and  its 
comparison  with  the  other  common  refrigerants  are  summarized 
in  Appendix  A,  pressure-enthalpy  diagram  of  R-114  is 
provided  in  Appendix  B  and  the  range  of  application  of 
refrigerants  are  provided  in  Appendix  C  [Ref.  4]. 

C.   ENHANCED  BOILING  TUBE  SURFACES 

The  evalution  of  enhanced  surface  geometries  that 
promote  high-performance  nucleate  boiling  has  been  surveyed 
by  Webb  [Ref.  5].  Two  basic  types  of  enhanced  boiling 
surfaces  are  employed  on  the  outer  surface  of  tubes:  porous 
coatings  and  reentrant  grooves.  The  High-Flux  surface  is  a 
sintered,  porous  metallic  matrix  bonded  to  the  base  surface. 
The  Gewa-T  and  Thermoexcel-E  surfaces  have  reentrant- type 
nucleation  sites  formed  by  cold  working  an  integral  finned 
tube  wall.  See  Figures  1.1,  1.2  and  1.3  for  details. 

The  key  to  the  high  performance  of  these  surfaces  is 
attributed  to  three  factors  [Ref.  5]  :  "(1)  a  pore  or  reent- 
rant cavity  within  a  critical  size  range,  (2)  interconnected 
cavities,  and  (3)  nucleation  sites  of  a  reentrant  shape. 
When  the  cavities  are  interconnected,  one  active  cavity  can 
activate  adjacent  cavities.  It  appears  that  the  dominant 
fraction  of  the  vaporization  occurs  at  a  very  thin  liquid 
film  within  the  subsurface  structure.  The  reentrant  cavity 
shape  provides  a  stable  vapor  trap,  which  will  remain  active 
at  very  low  liquid  superheat  values . " 


15 


The  pool  boiling  heat  -  transfer  performance  of  R-113  from 
the  above-mentioned  commercially-available  enhanced  surfaces 
has  been  investigated  by  Marto  and  Lepere  [Ref.  6].  They 
reported  that  for  R-113,  the  enhanced  surfaces  showed  a  two 
to  tenfold  increase  in  the  heat- transfer  coefficient  when 
compared  to  a  plain  tube.  The  High-Flux  surface  was  most 
effective  over  a  broad  range  of  heat  fluxes,  whereas  the 
Gewa-T  surface  was  only  effective  at  high  heat  fluxes  (near 
100  kW/m2 ) .  The  Thermoexcel-E  surface  showed  similar  gains 
in  the  heat- transfer  coefficient  to  that  of  the  High-Flux 
surface  below  10  kW/m2 .  They  also  reported  that  the  degree 
of  superheat  required  to  activate  the  enhanced  surfaces  was 
less  than  the  plain  tube,  and  was  sensitive  to  initial 
surface  conditioning  and  fluid  properties. 

D.   NUCLEATE  POOL  BOILING  PERFORMANCE  OF  PURE  R-114 

The  heat- transfer  performance  of  pure  R-114  in  .the 
evaporator  section  of  an  air-conditioning  unit  can  be  esti- 
mated using  one  of  the  suggested  commonly-.used  correlations 
for  nucleate  pool  boiling,  which  are  listed  by 
Chongrungreong  and  Sauer  [Ref.  7].  However,  because  of  the 
lack  of  complete  understanding  of  the  nucleate  boiling 
mechanism,  these  correlations  should  be  compared  with  exper- 
imental data  if  it  is  possible.  So  far,  the  literature 
reveals  a  limited  amount  of  nucleate  boiling  experimental 
investigations  using  R-114.  Happel  [Ref.  8],  using  a  hori- 
zontal, pure-nickel  tube  with  a  surface  roughness  of  0.43 
micrometer,  investigated  the  boiling  characteristics  of  pure 
R-114  at  various  pressures.  The  effect  of  pressure  on 
nucleate  boiling  heat  transfer  in  R-114  was  also  investi- 
gated experimentally  by  Nishikawa  and  his  co-workers 
[Ref.  9].  They  operated  over  a  reduced  pressure  (P/Pc  ) 
range  from   0.03  to   0.98,   and   demonstrated  that   the  heat 
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transfer  coefficient  of  R-114  increases  remarkably  at  high 
pressure  near  the  critical  point.  Their  results  showed  that 
the  effect  of  surface  roughness  changes  with  pressure  and 
the  effect  gradually  vanishes  as  pressure  approaches  the 
critical  point.  Using  a  horizontal  smooth  copper  tube, 
Henrici  [Ref.  10]  investigated  boiling  heat- transfer 
performance  of  R-114  over  a  heat-flux  range  from  0.1  to  100 
kW/m2  and  over  a  pressure  range  from  37  to  252  kN/m2 . 

Stephan  and  Mitrovic  [Ref.  11]  investigated  the  perform- 
ance of  pure  R-114  from  a  Gewa-T  enhanced  surface.  Boiling 
occurred  outside  a  horizontal  bundle  of  these  tubes.  They 
reported  that  the  position  of  the  tube  within  the  bundle 
does  not  effect  the  heat- transfer  coefficient  considerably. 
The  reason  for  this  effect  was  explained  in  the  following 
manner:  "The  flow  pattern  around  a  tube  seems  to  be  of  very 
weak  influence  and  the  heat- transfer  process  is  mainly 
governed  by  the  bubble  formation  and  the  two-phase  flow 
inside  the  T-shaped  channel.  The  top  part  of  the  fins,  which 
were  rolled  to  form  the  T,  practically  did  not  directly 
contribute  to  the  heat  -  transfer  process  because,  even  at 
maximum  heat  fluxes,  vapour  bubbles  were  not  produced 
there. " 

Marto  and  Hernandez  [Ref.  12]  also  investigated  nucleate 
pool-boiling  characteristics  of  a  Gewa-T  surface  in  R-113. 
Their  experimental  results  showed  a  disagreement  with 
Stephan  and  Mitrovic 's  above  explanation.  They  showed  that 
"the  liquid-vapor  motion  may  not  always  be  confined  to  flow 
circumf erentially  within  the  channels  formed  by  the  T-shaped 
fins  as  previously  postulated.  Surface  imperfections  on  the 
tips  o-f  the  Gewa-T  fins  may  cause  bubble  nucleation  from 
these  sites.  The  presence  of  these  active  sites  may  further 
enhance  the  performance  of  the  Gewa-T  surface." 

The  factors  which  appear  to  affect  boiling  heat- transfer 
performance  of  a  single  tube   are  the  physical  properties  of 
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refrigerant,  boiling  pressure,  tube  diameter,  surface  condi- 
tion of  the  tube,  and  maybe  the  "submergence  effect" 
(resulting  from  the  hydrostatic  liquid  head)  on  the  evapora- 
tion process.  Experimental  observations  have  shown  that 
changes  in  magnitude  of  these  properties  and  conditions 
significantly  affect  pool-boiling  heat-transfer  performance. 
[Ref.  7]. 

E.   NUCLEATE  BOILING  PERFORMANCE  OF  REFRIGERANT-OIL  MIXTURES 

Most  of  the  refrigeration  systems  work  with  oil- 
lubricated,  hermetically  sealed  compressors,  which  are 
cooled  by  the  working  fluid.  Therefore,  the  fluids  being 
circulated  within  these  refrigeration  systems  generally  are 
not  pure  refrigerants,  but  are  refrigerant-oil  mixtures. 
The  saturation  temperature  of  a  ref rigerant- oil  mixture  is 
higher  than  that  for  the  pure  refrigerant  at  the  same  pres- 
sure. Also,  both  viscosity  and  surface  tension  increase 
when  oil  is  added  to  refrigerants.  Other  thermodynamic  prop- 
erties of  refrigerants,  such  as  density  and  specific  heat, 
are  also  significantly  changed  by  the  addition  of  oil.  Thus, 
adequate  knowledge  on  the  influence  of  oil  on  the  boiling  of 
refrigerants  is  very  important  for  the  successful  design  of 
refrigeration  systems.  The  literature  reveals  only  very 
little  information  on  the  boiling  performance  of  R-114-oil 
mixtures  [Ref.  10].  In  particular,  heat- transfer  perform- 
ance of  the  enhanced  evaporator  surfaces  in  refrigerant-oil 
mixtures  is  a  poorly-understood  subject  area.  Since  these 
advanced  boiling  surfaces  contain  a  high  density  of  pores, 
the  presence  of  oil  can  generally  result  in  a  greater 
performance  degradation  than  on  more- conventional ,  smooth 
surfaces . 

Various  studies  which  are  surveyed  and  listed  by  Jensen 
and  Jackman  [Ref.  13]  have  shown  that,  in  general,   nucleate 
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boiling  heat- transfer  coefficients  decrease  with  increasing 
oil  concentration  in  R-ll,  R- 12 ,  R-22,  R-113  and  R-114.  For 
oil  concentrations  less  than  3%,  however,  Dougherty  and 
Sauer  [Ref.  14]  demonstrated  that  the  heat-tansfer  coeffi- 
cient may  be  increased  slightly  for  R-ll-oil  and  R-113-oil 
mixtures.  They  attributed  this  trend  to  an  increase  of 
foaming  of  the  mixtures.  This  result  also  has  been • confirmed 
by  Sauer,  Gibson  and  Chongrungreong  [Ref.  15]  using  R-12-oil 
mixtures,  and  by  Henrici  [Ref.  10]  using  R-114-oil  mixtures. 
However,  in  a  recent  experimental  study  using  R-113-oil 
mixtures,  Jensen  and  Jackman  [Ref.  13]  did  not  observe  such 
an  enhancement  with  low  oil  concentrations.  The  model  of 
the  governing  processes  involved  in  boiling  of  refrigerant- 
oil  mixtures  has  been  given  by  Stephan  [Ref.  16]  as  follows: 
"Obviously  the  reduction  in  heat  transfer  in  a  mixture  is 
due  to  the  fact  that  the  liquid  near  the  vapour- liquid 
interface  of  growing  bubbles  becomes  denuded  in  the  more 
volatile  components.  Thus,  a  concentration  gradient  exists 
between  the  interface  and  the  bulk  of  the  liquid,  and  the 
more  volatile  components  have  to  diffuse  through  the  liquid 
in  order  to  reach  the  interface.  A  mass- transfer  resistance 
thus  is  introduced,  which  is  not  present  when  boiling  single 
components . " 

The  effect  of  the  oil  viscosity  on  heat  transfer  in 
boiling  of  refrigerant-oil  mixtures  has  been  considered  by 
Sauer,  Gibson,  and  Chongrungreong  [Ref.  15].  Oil  viscosity 
is  believed  to  significantly  contribute  to  the  changes  in 
boiling  characteristics  of  ref rigerant -oil  mixtures.  The 
heat- transfer  coefficient  is  greatest  for  the  highest 
viscosity  of  oil  for  the  same  temperature  difference. 

The  effect  of  oil  on  boiling  of  R-12  from  a  Gewa-T  tube 
was  investigated  by  Stephan  and  Mitrovic  [Ref.  17].  They 
reported  that  the  presence  of  oil  in  R-12  reduced  the  heat- 
transfer  coefficient   almost  over  the   entire  range   of  heat 
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fluxes  and  mass  fractions .  They  observed  an  increase  in 
heat  -  trans fer  coefficient  with  less  than  5%  mass  fraction  of 
oil  only  at  high  heat  fluxes  (such  as  higher  than  20  kW/m2 ) . 
The  effects  of  oil  on  the  heat  -  transfer  performance  of 
refrigerants  in  a  tube  bundle  have  been  summarized  by  Arai , 
et  al .  [Ref.  18]  in  the  following  manner:  "The  presence  of 
oil  in  a  shell-and- tube  evaporator  is  believed  to  have  two 
opposing  effects  on  the  heat- transfer  coefficient  of  the 
evaporator  tube  bundles.  Firstly,  many  people  suspect  that 
the  oil  left  out  from  the  evaporating  refrigerant  tends  to 
choke  the  surface  structure,  thereby  lowering  the  heat- 
transfer  coefficient.  Secondly,  from  the  experience  of 
conventional  heat  exchangers,  it  has  long  been  recognized 
that  the  presence  of  oil  enhances  foaming  on  the  upper  part 
of  the  tube  bundles.  This  foaming  increases  the  vigor  of 
refrigerant  movement  around  the  tubes,  so  that  heat  transfer 
is  enhanced . " 

F.   THESIS  OBJECTIVES 

Based  on  the  foregoing  discussion,  it  is  clear  that 
little  information  exists  in  the  literature  to  enable  the 
successful  design  of  R-114  evaporators,  especially  when 
using  enhanced  boiling  surfaces,  in  the  presence  of  lubri- 
cating oils.  The  objectives  of  this  research  effort  were 
therefore  to: 

1.  Design,  construct  and  instrument  an  experimental 
apparatus  in  order  to  investigate  the  heat- transfer 
performance  of  R-114  and  R-114-oil  mixtures  from 
enhanced  evaporator  tube  surfaces , 

2.  Successfully  operate  the  system  to  produce  repeatable 
data,  and 

3.  Obtain  smooth-tube  baseline  data  for  R-114  to  be  used 
as  a  standard  for  the  future  data  on  advanced 
surfaces . 


20 


GEWA-Tube 


iigsiiissp 


.?'■•■''     if?       f      I   #-» 


Figure  1.1    Schematic  and  Optical  Micrograph  of  a  Cross 
Section  of  the  Gewa-T  Surface  I20X). 
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Figure    1 . 2 


Schematic  and  Optical  Micrograph  of  the 
Thermoexcel-E  Surface  (40X). 
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Figure    1 . 3 


Schematic  and  Scanning  Electron  Micrograph  of  the 
High  Flux  Surface  (500X). 
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II.  DESCRIPTION  OF  EXPERIMENTAL  APPARATUS 

A.   OVERVIEW  OF  THE  SYSTEM 

An  overall  schematic  representation  of  the  experimental 
apparatus  is  shown  in  Figure  2.1,  and  a  photograph  is  shown 
in  Figure  2.2.  The  apparatus  consists  of  a  Pyrex-glass 
evaporator  with  the  boiling  test  tube  1,  a  Pyrex-glass 
condenser  2,  an  R-114  liquid  reservoir  3,  an  oil  reservoir 
4,  a  graduated  oil  cylinder  5,  a  vacuum  pump  6,  a  cooling 
section  consisting  of  a  1/2-Ton  R- 12  refrigeration  unit,  a 
tank  filled  with  water-ethy lene  glycol  mixture  7,  and  a  pump 
8.  The  evaporator  1,  condenser  2,  oil-reservoir  4, 
graduated-oil  cylinder  5  and  R-114  reservoir  3,  are  placed 
in  a  Plexiglas  chamber  9,  where  the  temperature  can  be  main- 
tained at  lower  values  (up  to  10  °C)  relative  to  ambient 
temperature.  This  chamber  also  served  as  a  safety  barrier  in 
case  of  an  inadvertent  overpressurization  of  the  glass 
apparatus . 

The  vapor  rising  from  the  evaporator  is  condensed  in  the 
condenser  2  and  the  liquid  is  then  fed  back  into  the  evapo- 
rator through  valve  V5 .  In  this  manner,  the  specific 
enthalpy  of  condensation  ("latent  heat")  of  R-114  vapor  is 
removed  in  the  condenser.  The  evaporator  of  the  R-12 
refrigeration  system  is  made  of  a  copper  coil,  which  is 
immersed  in  the  water-ethylene  glycol  tank  in  order  to  main- 
tain the  temperature  of  the  mixture  at  about  -17  °C.  An 
8-GPM,  turbine-type  pump  8  transports  the  water-ethylene 
glycol  mixture  through  the  control  valve  VC  to  the 
condenser. 

The  oil  can  be  added  into  the  evaporator  1  using  the 
graduated   cylinder  5   by   opening  valves   VI   and  V4 .    The 
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graduated  cylinder  is  supplied  with  oil  by  reservoir  4.  The 
valves  V5 ,  V6  and  V7  are  provided  for  emptying  the  evapo- 
rator prior  to  dismantling  it,  for  example,  for  changing  the 
evaporator  tube.  Also,  fresh  R-114  can  be  added  to  the 
reservoir  while  valve  V5  is  closed  and  V6 ,  V7  are  open  and 
by  providing  R-114  vapor  through  V8  and  V10 . 

The  test  tube  was  placed  horizontally  in  the  evaporator 
after  a  240  VAC  electrical  cartridge  heater  of  203.2  mm  (8 
in.)  in  length  had  been  inserted  in  its  middle  section.  The 
whole  system  was  connected  to  the  vacuum  pump  through  V8  and 
V10  in  order  to  remove  noncondensible  gases.  The  relief 
valve  was  set  to  a  gage  pressure  limit  of  138  kN/m2  (20 
psi),  which  is  50%  less  than  the  manufacturer-recommended 
working  pressure  limit  of  the  Pyrex  glass  tees.  All  connec- 
tions were  assembled  with  Swagelok  fittings  sealed  with 
Teflon  ferrules  to  ensure  leak  tightness.  Also,  all  valves 
used  in  the  system,  except  control  valve  VC  and  by-pass 
valve  V9 ,  were  made  of  stainless  steel.  Copper  tubes  of 
9.5-mm  (3/8  in.)  diameter  were  used  for  all  piping,  except 
in  the  oil  section,  where  6.35-mm  (1/4  in.)  tubes  were  used. 

B.   BOILING  TEST  SECTION 

1 .   Evaporator 

The  evaporator  is  a  T-shaped  container  made  of 
Corning  Pyrex  glass.  Figure  2.3  shows  a  schematic  of  the 
evaporator  with  dimensions.  The  selection  of  Pyrex  glass 
offers  several  advantages:  it  is  corrosion- resistant  and 
transparent  and  has  a  smooth  interior  surface  (this  mini- 
mized nucleate  boiling  at  the  inner  surface  of  the 
container)  compared  to  metals,  and  is  more  stable  with 
temperature  and  pressure  variations  compared  to  ordinary 
glass.  The  operating  pressure  limit  of  the  evaporator  is 
guaranteed  by  the  manufacturer  up  to   a  gage  pressure  of  200 
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kN/m2  (30  psi).  Each  end  of  the  Pyrex-glass  evaporator  was 
fitted  with  a  cast-iron  flange  and  a  gasket.  A  detailed 
sketch  of  the  cast-iron  flange  is  shown  in  Figure  2.3.  Two 
aluminum  flanges,  210  mm  in  diameter  and  12.7  mm  in  thick- 
ness, were  bolted  to  the  cast-iron  flanges.  Thus,  all 
fittings  were  connected  through  the  aluminum  flanges  to  the 
Pyrex-glass  evaporator. 

2 .  Test  Section 

A  schematic  drawing  of  the  test  section  is  shown  in 
Figure  2.4.  The  boiling  test  tube  was  held  in  place  by  two 
Teflon  bushings,  which  were  attached  to  the  aluminum  flanges 
at  both  ends  of  the  T-shaped  evaporator.  Four  studs  were 
used  in  order  to  attach  the  Teflon  bushing  on  the  aluminum 
flange.  The  Teflon  bushing  and  both  ends  of  the  test  tube 
were  sealed  by  means  of  two  rubber  0-rings  in  order  to 
achieve  leak  tightness.  An  additional  0-ring  was  placed 
between  the  aluminum  flange  and  the  Teflon  plug  at  each  end. 

3 .  Boiling  Tube 

A  smooth,  hard-copper  tube,  15.9  mm  (5/8  in.)  in 
outer  diameter,  12.7  mm  (1/2  in.)  in  inside  diameter  and 
431.8  mm  (17  in.)  in  length,  was  used  to  provide  baseline 
data  for  the  nucleate  pool  boiling  heat  -  transfer  coefficient 
of  pure  R-114.  The  heater  was  a  1000-Watt  240-Volt 
stainless- steel  cartridge,  6.35  mm  (1/4  in.)  in  outer  diam- 
eter and  203.2  mm  (8  in.)  in  length.  The  cartridge  heater 
was  inserted  into  a  copper  sleeve,  which  was  6.35  mm  (1/4 
in.)  in  inside  diameter,  12.7  mm  (1/2  in.)  in  outside  diam- 
eter and  203.2  mm  (8  in.)  in  length.  In  order  to  provide  a 
uniform  heat  flux,  the  cartridge  heater  and  the  copper 
sleeve  were  soldered  together.  Pulido  [Ref.  19]  showed  that 
(using  a  similar  type  of  heater-element  construction)  the 
heater   element   provides   a  uniform   heat   flux   along   the 
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circumferential  direction.  The  heater  element  used  during 
this  thesis  was  assumed  to  have  similar  characteristics. 
The  cartridge  heater  and  the  copper  sleeve  were  then 
inserted  as  a  unit  into  the  middle  portion  of  the  test  tube 
using  a  tight  mechanical  fit.  (Note:  A  shrink-fit  process 
is  essential  to  minimize  the  thermal  contact  resistance. 
However,  owing  to  the  time  constraint  in  this  thesis,  the 
shrink-fit  process  was  delayed  and  a  slide  fit  with  a  clear- 
ance of  less  than  0.01  mm  was  used.  Thus,  the  resulting  tube 
was  used  mainly  for  the  purpose  of  obtaining  a  successfully- 
operating  apparatus).  The  active  boiling  length  of  the  test 
tube  was  203.2  mm  (8  in.)  in  the  middle  portion  of  the 
boiling  tube.  Tn  order  to  compute  the  actual  average  heat 
flux  in  the  heated  portion  of  the  tube,  a  suitable  correc- 
tion was  applied  for  the  natural  convection  created  at  both 
ends.  See  Heat  Flux  Calculation  section  in  Chapter  4  for 
details . 

To  measure  the  wall  surface  temperature  of  the 
boiling  tube,  8  thermocouples  were  inserted  into  8  grooves 
which  were  machined  on  the  outside  of  the  copper  sleeve.  As 
shown  in  Figure  2.5,  these  thermocouples  were  located  at 
different  axial  and  circumferential  locations.  Four  thermo- 
couples (1,2,5,6)  were  placed  with  90-degree  separation 
around  the  copper  sleeve  to  measure  the  temperature  distri- 
bution at  the  mid  cross  section  of  the  active  length  of 
boiling  tube.  The  thermocouple's  wire  was  glued  down  at 
several  places  along  the  thermocouple's  channel  using  Epoxy. 
The  longitudinal  temperature  distribution  on  the  active 
boiling  section  was  also  measured  by  the  symmetrically- 
located  4  thermocouples  (3,4,7,8).  The  exact  locations  of 
the  thermocouples  and  dimensions  of  the  thermocouple  grooves 
are  given  in  Figure  2.5.  All  thermocouple  grooves  were 
axially  machined  from  the  location  of  the  thermocouple  hot 
junctions  to  the  nearest  end  of  the  copper  sleeve. 
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C.  CONDENSER  SECTION 

The  condenser  was  also  a  T-shaped  container  made  of 
Corning  Pyrex  glass.  It  was  identical  to  the  evaporator.  The 
position  of  the  condenser  can  be  seen  in  Figure  2.1.  R-114 
was  condensed  on  a  vertical  copper  coil,  which  was  inserted 
in  the  Pyrex-glass  condenser.  The  copper  coil  was  fabricated 
into  a  76.2  mm  (3  in.)  diameter  coil  using  9.5  mm  (3/8  in.) 
copper  tube.  The  active  condensation  length  was  estimated  to 
be  4.5  m  (15  ft) . 

The  top  portion  of  the  condenser  was  connected  to  a 
portable,  mechanical  vacuum  pump  to  remove  noncondensible 
gases  from  the  apparatus.  The  bottom  of  the  condenser  was 
also  connected  to  the  evaporator  via  valve  V5  in  order  to 
return  the  condensed  R-114  liquid  to  the  evaporator.  The 
cooling  liquid,  i.e.,  water-ethylene  glycol  mixture,  entered 
the  top  portion  of  the  condenser,  through  the  copper  coil 
and  left  the  condenser  from  the  bottom  to  return  to  the 
water-ethylene  glycol  tank.  The  condenser  was  placed  verti- 
cally and  connected  to  the  vapor  outlet  of  the  evaporator 
using  L-shaped  aliminum  tube,  50.8  mm  (2  in.)  in  diameter. 
The  maximum  vapor  velocity  of  R-114  vapor  through  the 
aluminum  tube  was  found  to  be  about  0.5  m/sec.  A  Bourdon 
gage  with  a  range  of  absolute  vacuum  to  a  gage  pressure  of 
1030  kN/m2  (150  psi)  and  a  relief  valve  which  was  set  to  138 
kN/m2  (20  psi)  were  placed  on  the  L-shaped  aluminum  tube. 

D.  OIL  ADDING  SECTION 

To  study  the  boiling  performance  of  R-114-oil  mixtures, 
a  cylindrical  aluminum  reservoir,  152.4  mm  (6. in)  in  diam- 
eter and  152.4  cm  (6  in.)  in  height,  and  a  glass  oil 
cylinder  were  installed  above  the  evaporator.  The  relative 
positions  of  the  oil  reservoir  and  the  oil  cylinder  are 
shown  in  Figure  2.1.    The  oil  cylinder  was  355  mm  in  length 
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and  had  a  diameter  of  25.4  mm.  This  cylinder  was  specially 
ordered  to  achieve  a  resolution  of  0.5  ml.  The  oil  cylinder 
was  connected  to  the  oil  reservoir  through  valves  V3  and  V2 . 
The  addition  of  oil  into  the  evaporator  can  be  achieved 
through  VI  by  gravity  after  balancing  the  pressure  of  the 
oil  cylinder  with  that  in  the  evaporator  by  opening  valve 
V4. 

E.   COOLING  SECTION 

1.  Water-Ethylene  Glycol  Mixture  Tank 

In  order  to  store  the  water-ethy lene  glycol  mixture, 
a  special  tank  was  manufactured.  The  total  volume  of  the 
tank  was  0.154  m3  (0.48  m  x  0.48  m  x  0.66  m)  and  it  was  made 
of  12 . 7-mm- thick  Plexiglas  sheet.  All  sides  of  the  tank  were 
glued  together  with  methylene- chloride  solution.  The  joints 
were  held  together  with  small  screws  for  extra  strength.  The 
low  thermal  conductivity  of  Plexiglas  was  especially  suited 
to  minimize  heat  transfer  (from  room  to  water-ethylene 
glycol  mixture)  through  the  tank  walls.  The  tank  was  placed 
on  the  floor  and  all  sides  were  insulated  with  22  mm  (7/8 
in.)  thick  rubber  insulation  sheets.  The  cooling  mixture 
contained  49  liters  (13  gal.)  of  ethylene  glycol  and  94 
liters  of  (25  gal.)  distilled  water.  The  freezing  point  of 
this  mixture  was  about  -  25  °C. 

2 .  R- 12  Refrigeration  Plant 

A  1/2-Ton  R- 12  refrigeration  plant  was  installed  to 
cool  the  water-ethylene  glycol  mixture.  Figure  2.6  shows  a 
schematic  of  the  R- 12  refrigeration  plant.  It  consists  of  a 
compact-type  air-cooled  condenser,  a  compressor,  a  receiver, 
a  filter-drier  unit,  a  pressure  regulator,  a  pressure- 
control  switch  and  a  thermostatic  expansion  valve.  The 
evaporator  of   the  R- 12  refrigeration  plant   was  constructed 
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using  a  9.5  mm  (3/8  in.)  copper  tube,  which  was  immersed  in 
the  water- ethylene  glycol  tank.  The  temperature  of  the 
water- ethylene  glycol  mixture  is  controlled  by  both  a  ther- 
mostatic expansion  valve  and  a  pressure  control  switch.  The 
R- 12  refrigeration  plant  was  adjusted  to  keep  the  tempera- 
ture of  the  cooling  liquid  at  about  -17  °C. 

3 .   Pump  and  Control  Valve 

An  8  GPM,  115  VAC  Burks  turbine- type ,  positive- 
displacement  pump  was  installed  on  the  floor  and  the 
suction  side  of  the  pump,  25.4  mm  (1  in.)  in  diameter,  was 
directly  coupled  to  the  water-ethy lene  glycol  tank.  Cooling 
liquid  was  pumped  from  the  tank  to  the  condenser  through  the 
control  valve  VC .  Also,  a  by-pass  valve  V9  was  placed  before 
the  control  valve  VC  on  the  discharge  line.  The  use  of  the 
by-pass  line  served  two  important  purposes:  1,  it  avoided 
overloading  of  pump  8  in  the  event  value  VC  is  completely 
closed,  and  2,  it  provided  proper  mixing  for  the  "warm" 
stream  returning  from  the  condenser  2.  See  Figure  2.1  for 
the  positions  of  control  valve  VC  and  by-pass  valve  V9 .  The 
by-pass  valve  was  adjusted  and  set  permanently  so  that  a 
sufficient  range  of  flow  rates  can  be  achieved  through  valve 
VC  to  satisfy  proper  cooling  of  the  condenser. 

F.   R-114  RESERVOIR 

An  aluminum  cylindrical  reservoir,  228  mm  (9  in.)  in 
diameter  and  254  mm  (10  in.)  in  height,  was  placed  verti- 
cally between  the  evaporator  and  condenser  in  order  to  store 
R-114  as  a  liquid.  The  liquid  level  of  the  R-114  can  be 
observed  by  means  of  a  sight  hose  attached  on  the  reservoir 
with  proper  fittings.  R-114  reservoir  was  connected  to  the 
vapor  line  through  valve  V7  and  to  the  liquid  line  through 
valve  V6 .  See  Figure  2.1  for  arrangement  of  the  reservoir. 
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G.   CHAMBER 

An  aluminum  frame  (1.07  mx  0.51m  x  0.61  m)  was 
constructed  to  locate  all  the  parts  of  the  apparatus,  except 
the  cooling  section.  All  four  vertical  sides  of  the  frame 
were  covered  with  12.7  mm  (1/2  in.)  thick  Plexiglas  sheets 
and  both  left  and  right  sides  were  provided  with  hinges  to 
enable  easy  access  to  the  components  of  the  apparatus. 
Aluminum  and  plywood  plates  were  used  to  cover  the  bottom 
and  top  sides  of  the  frame,  respectively.  The  valve  bodies 
of  VI  through  V8  were  placed  inside  of  the  front  Plexiglas 
sheet  with  the  valve  stems  penetrating  the  sheet.  Thus,  the 
valve  handles  were  accessible  from  outside  of  the  Plexiglas 
box.  The  whole  frame  was  placed  above  the  water- ethylene 
glycol  tank  with  aluminum  support  so  that  the  system  was 
very  compact . 

One  of  the  main  advantages  of  this  chamber  is  that  the 
temperature  surrounding  the  evaporator  can  be  reduced  rela- 
tive to  the  ambient  temperature.  Also,  in  case  of  emergency, 
the  thick  Plexiglas  chamber  would  provide  a  safety  barrier 
to  personnel  and  equipment . 

H.   INSTRUMENTATION 

1 .   Power  Measurement 

A  240  Volt  AC  source  was  used  as  the  power  supply, 
and  it  was  adjusted  by  a  variac  in  the  range  of  0-260  Volt 
and  0-8  Ampere  according  to  the  desired  heat  flux  at  the 
surface  of  the  boiling  tube.  Power  input  to  the  boiling  tube 
was  measured  with  an  AC  current  sensor  and  an  AC-DC  true 
R.M.S  converter.  The  AC  current  sensor  was  connected  to  the 
input  line  of  the  heater  in  series  and  the  AC-DC  true  R.M.S 
converter  was  connected  in  parallel.  Figure  2.7  shows  a 
schematic  representation   of  the   power-measurement  devices. 
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Both  the  AC  current  sensor  and  the  AC-DC  true  R.M.S 
converter  were  connected  to  the  data  acquisition/ reduction 
system . 

2 .   Temperature  Measurement 

Various  temperatures  were  monitored  throughout  the 
system  to  include: 

1.  Boiling  tube  wall  (8  thermocouples) 

2.  Liquid  temperature  (one  thermocouple) 

3.  Vapor  temperature  (two  thermocouples),  and 

4.  Water- ethylene  glycol  mixture   temperature  (one  ther- 
mocouple ) 

The  locations  of  the  wall  thermocouples  in  the  sleeve  of  the 
boiling  tube  are  shown  in  Figure  2.5.  The  liquid  and  vapor 
thermocouples  were  inserted  into  the  two  specially- 
manufactured  thermocouple  wells.  Figure  2.8  shows  a  sche- 
matic of  these  thermocouple  wells.  While  the  stainless- steel 
portion  minimizes  (owing  to  low  thermal  conductivity)  errors 
resulting  from  the  axial  conduction  of  heat  from  the 
surrounding,  the  copper  tip  helps  minimize  the  temperature 
drop  from  the  area  being  measured  to  the  thermocouple  loca- 
tion (owing  to  the  high  thermal  conductivity  of  copper). 

All  the  temperature  measurements  were  accomplished 
by  0.245  mm  (30  gage)  copper- constantan  thermocouples.  Each 
thermocouple  measurement  was  read  directly  by  a 
Hewlett-Packard  3497A  data  acquisition  system,  which  was 
controlled  by  a  Hewlett-Packard  9826  computer.  Each  thermo- 
couple was  scanned  for  0.8  seconds  and  twenty  readings  were 
averaged  to  obtain  a  more  accurate  measurement. 

A  total  of  five  thermocouples  were  calibrated.  Two 
thermocouples  were  made  from  the  beginning  of  a  spool  of 
copper- constantan  wire;  one  was  made  from  the  mid  portion 
and  two  were  from  the  end.  It  was  assumed  that  the  proper- 
ties of   the  copper   and  constantan   do  not   change  along   a 
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given  section  of  wire  for  any  given  spool.  All  thermocouples 
were  calibrated  by  the  method   described  in  Appendix  D. 
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Figure  2.1    Schematic  of  the  Boiling  Test  Apparatus 
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Figure  2.2    Photograph  of  Overall  Syste 
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D, 


a)  Corning  Pyrex  Glass  Evaporator  (D  x  Di 
L  =  17  8  mm ,  L]_  =  12  7  mm ) 


4  0  2x51  mm, 


L.T^^TTT 


3k. 


INSERT 


f 


b)  Cast  Iron  Flange  and  Gasket  (dj=  190mm,  d0 =  210mm, 
L  =  14mm,  A  =  2 1°  ) 

Figure  2.3    Schematic  of  the  Pyrex  Glass  Evaporator. 
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III.  DATA  ACQUISITION/ REDUCTION 

A.  DATA  ACQUISITION  AND  STORAGE 

A  Hewlett-Packard  3497A  automatic  data  acquisition/ 
control  unit  was  used  to  read  temperatures  from  the  thermo- 
couples and  to  read  current  and  voltage  values  of  the 
cartridge  heater  from  the  AC  current  sensor  and  the  AC-DC 
true  RMS  converter,  respectively.  A  Hewlett-Packard  9826A 
computer  was  used  to  control  the  data  acquisition/ control 
unit  and  to  analyze  and  store  data.  Figure  3.1  shows  a 
photograph  of  the  data  acquisition/ control  unit  and 
computer. 

Information  was  entered  through  the  keyboard  to  prompt 
the  data  acquisition/ control  unit  to  automatically  scan  each 
channel.  Channel  assignments  are  listed  in  Table  I.  These 
raw  data  were  immediately  processed  and  a  hard-copy  printout 
was  provided.  Also,  these  data  were  transferred  to  a 
computer  disk  under  a  user- specif ied  file  name  to  keep  a 
permanent  record.  The  ability  to  store  raw  data  directly 
enabled  these  data  to  be  reduced  at  any  time  and  allowed 
flexibility  for  changes  to  data- reduction  software. 

B.  DATA  REDUCTION 

Following  data  acquisition  for  each  data  point,  results 
were  computed  according  to  the  stepwise  procedure  outlined 
in  the  next  section,  and  then  printed  on  a  Hewlett-Packard 
2671G  thermal  printer.  Heat  flux  versus  wall  superheat  (the 
temperature  of  the  boiling  surface  minus  the  liquid  tempera- 
ture) were  also  stored  in  a  user- specif ied  plot  data  file 
for  subsequent  plotting  using  the  subprogram  PLOT.  A 
Hewlett-Packard  7470A  plotter  was  also  interfaced  with  the 
computer . 
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STEPWISE  DATA- COLLECTION  AND  SOLUTION  PROCEDURE 

1.  Select  tube  type  (all  dimensions  of  the  boiling  test 
tube  are  included) . 

2.  Set  desired  heat  flux  and  saturation  temperature  of 
the  boiling  liquid.  Wait  for  steady-state  conditions 
(See  Chapter  4  section  B  for  details). 

3.  Scan  all  channels  listed  in  Table  I  (thermocouples, 
AC  current  sensor  and  AC-DC  true  RMS  converter  read- 
ings) . 

4.  Save  these  raw  data  in  a  user- specif ied  file. 

5.  Convert  these  raw  data  readings  to  corresponding 
units  (temperature  in  °C,  current  in  Ampere,  voltage 
in  volt ) . 

6.  Compute  the  heat- transfer  rate  from  the  cartridge 
heater  (See  Chapter  4,  Section  C  for  details). 

7.  Compute  the  average  wall  temperature  of  the  boiling 
test  tube  and  calculate  the  wall  superheat  of  this 
data  set  (See  Chapter  4,  Section  C  for  details). 

8.  Compute  the  physical  properties  of  R-114  using  given 
correlations  at  film  temperature  (See  Appendix  H  for 
details ) . 

9.  Compute  the  natural- convection  heat- transfer  coeffi- 
cient of  R-114  from  non-boiling  ends  of  the  test 
tube  . 

10.  Compute  heat  losses  from  non-boiling  ends. 

11.  Calculate  the  heat  flux  from  boiling  test  tube  to  the 
boiling  liquid. 

12.  Calculate  boiling  heat- transfer  coefficient  of  the 
R-114  from  test  tube. 

13.  Store  the  heat  flux  versus  wall  superheat  values  for 
each  data  set  in  a  user- specif ied  plot  file. 

14.  Use  the  subprogram  PLOT  and  plot  the  data  run. 
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Appendix  E  shows  a  listing  of  computer  program  DRPR1  and 
subprogram  PLOT  and  Appendix  F  shows  an  example  of  represen- 
tative data  run. 

TABLE  I 
HP  3497A  Channel  Assignments 

Channel  Assignment 


25  -  32  Boiling  tube  wall  temperature 

33  Liquid  temperature 

34  -  35  Vapor  temperature 

36  Temperature  of  cooling  liquid 

62  AC-DC  true  RSM  converter 

63  AC  current  sensor 
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IV.  EXPERIMENTAL  PROCEDURE 

A.   PREPARATION 

1.  Pressure  Test  of  the  Apparatus  with  Air 

Upon  assembling  all  the  components,  including  the 
boiling  test  tube,  the  experimental  apparatus  was  pressur- 
ized with  air  up  to  100  kN/m2  (15  psi).  A  soap-bubble  test 
was  first  carried  out  to  locate  any  leaks  from  each  compo- 
nent of  the  apparatus.  All  detected  leaks  were  succesfully 
and  systematically  fixed  at  the  end  of  each  pressure  test. 

2 .  Pressure  Test  of  the  Apparatus  with  R- 114 

A  second  pressure  test  was  performed  using  R-114 
vapor.  The  apparatus  were  charged  with  R-114  vapor  in  the 
following  manner.  The  apparatus  was  evacuated  to  25  in.  Hg . 
by  a  portable  mechanical  vacuum  pump  through  valves  Vll  and 
V8  (See  Figure  2.1  for  the  configuration  of  valves).  Upon 
closing  valve  Vll,  a  vapor  outlet  of  the  R-114  supply 
cylinder  was  connected  to  the  apparatus  by  means  of  valve 
V10 .  Thus,  the  pressure  of  the  apparatus  was  increased  up  to 
the  saturation  pressure  of  R-114  at  that  ambient  temperature 
(ambient  temperature  was  21  °C  and  the  corresponding  satura- 
tion pressure  (gage)  was  86  kN/m2  (12.5  psi)).  An  Automatic 
Halogen  Leak  Detector,  TIF  5000,  was  used  to  detect  R-114 
leakage  from  the  apparatus.  The  sensitivity  of  this  detector 
is  1/2  oz.  per  year,  3  PPM  max.  A  few  small  leaks  were 
observed  and  successfully  isolated. 

3.  Charging  the  Apparatus  with  R- 114 

Upon  completing  the  pressure  tests,  the  evaporator 
of  the  apparatus  was  filled  with  R-114  liquid  up  to  a  marked 
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position   (20  mm   above   the  boiling   test   tube)   with   the 
following  procedure  (See  also  Figure  2.1). 

1.  The  temperature  of  the  water-ethy lene  glycol  tank  was 
reduced  to  about  -17  °C  by  the  R-12  refrigeration 
unit  . 

2.  The  apparatus  was  evacuated  to   about  29  in  Hg  vacuum 
, by  the  vacuum   pump  (the  valves  Vll  and   V8  were  open 

and  V10  was  closed) . 

3.  Upon  closing  Vll,  the  vapor  outlet  of  the  R-114 
supply  cylinder  was  connected  to  V10 . 

4.  The  cooling  section  pump  8,  delivered  the  cooling 
liquid  from  the  water-ethylene  glycol  tank  to  the 
condenser  through  the  control  valve  VC . 

5.  The  R-114  vapor  condensed  on  the  condenser  coils  and 
the  liquid  was  collected  in  the  evaporator  by 
gravity,  while  V5  was  open  and  V6  was  closed. 

6.  Pump  8  was  then  stopped  when  the  desired  liquid  level 
of  the  evaporator  was  achieved. 

7.  The  pressure  of  the  apparatus  was  then  allowed  to 
increase  up  to  the  saturation  pressure  corresponding 
to  the  ambient  temperature. 

B.   NORMAL  OPERATION 

The  following  procedure  was  established  to  obtain  the 
heat- transfer  coefficient  of  R-114  from  the  smooth  copper 
test  tube: 

1.  The  test  tube  was  immersed  in  the  pool  (20  mm  below 
the  liquid  level). 

2.  The  R-12  refrigeration  unit  was  operated  24  hours  in 
advance  in  order  to  reduce  the  temperature  of  the 
water-ethylene  glycol  tank  to  about  -17  °C. 

3.  The  data  acquisition/ control  unit,  computer,  power 
supply  and  cooling  section  pump  were  switched  on. 
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4.  Any  accumulated  noncondensible  gases  were  evacuated 
by  the  portable  mechanical  vacuum  pump  through  valve 
V8  and  Vll  (see  Figure  2.1). 

5.  The  iteractive  computer  program  DRPR1  was  operated 
and  the  desired  heat  flux  and  saturation  temperature 
were  given  to  the  program  as  reference  values. 

6.  The  cartridge-heater  voltage  was  then  adjusted  by 
variac  in  order  to  maintain  the  desired  heat  flux. 
The  desired  and  actual  heat  fluxes  were  compared 
continiously  by  the  computer  program  until  they 
agreed  to  within  2  percent . 

7.  The  actual  versus  desired  saturation  temperature  of 
the  liquid  R-114  was  also  monitored  by  computer 
program.  The  amount  of  cooling  liquid,  which  was 
being  circulated  through  the  copper  coil  in  the 
condenser,  was  regulated  by  control  valve  VC  in  order 
to  obtain  nearly  constant  saturation  temperature 
(±0.2  °C)  at  a  given  heat  flux. 

8.  For  all  the  consecutive  settings  during  both 
increasing  and  decreasing  heat  fluxes,  the  boiling 
was  allowed  to  stabilize  for  five  minutes  at  each 
power  setting  and  raw  data  (thermocouple  readings,  AC 
current  sensor  readings  and  AC-DC  true  RMS  converter 
readings)  were  recorded  in  a  user- specif ied  file. 

9.  At  a  given  power  setting  and  a  saturation  tempera- 
ture, the  following  processed  data  were  recorded  as  a 
printout:  wall  temperatures  of  the  boiling  section, 
liquid  bulk  temperature,  vapor  temperature,  net  heat 
flux,  temperature  of  the  water- ethylene  glycol 
mixture,  wall  superheat  and  heat- transfer  coefficient 
of  the  R-114. 

10.  For  each  data  set,  the  above  procedure  beginning  with 
step  5  was  repeated. 
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11.  For  each  data  series,  seven  different  heat  fluxes  (1, 
2,  5,  10,  20,  50  and  95  kW/m2 )  were  selected. 

C.   HEAT -FLUX  CALCULATION 

According   to   the   description    of   the   boiling  tube 

construction,  the  cartridge  heater,   which  was  inserted  into 

the  boiling  tube  is  the  heat  source.   The  variations  in  heat 

flux   are  made   through  different   voltage   settings  of  the 
variac . 

(^  =   VI  (1) 


where 


QH  =   heat- transfer  rate  from  the  cartridge 
heater  (W) 

V   =   Voltage  across  the  cartridge  heater 
element  (volts) 

I   =   current  through  the  heater  element  (amps) 


This  circuit  is  connected  in  series  with  an  AC  current 
sensor  and  in  parallel  with  an  AC-DC  true  RMS  converter. 
Each  of  these  two  sensors  produces  a  DC  output  in  the  range 
of  0-10  V,  which  is  automatically  read  by  the  data 
acquisition/ control  unit.  The  calibration  equations  of  these 
sensors  give  the  actual  values  as  follows: 

V   =   25VS  (2) 

I   =   Is  (3) 
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where : 


V   =   voltage  reading  by  AC-DC  true  RMS 
converter  (volts) 


I   =   current  reading  by  AC  current 
sensor  (volts) 

The  total  length  of  the  boiling  test  tube  was  355.6  mm, 
while  the  active  boiling  length;  i.e.,  the  length  of  the 
cartridge  heater,  was  203.2  mm.  The  geometry  o-f  the  boiling 
test  tube  is  shown  schematically  in  Figure  4.1. 
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Figure  4.1    Geometry  of  the  Boiling  Test  Tube. 

Natural- convection  heat  transfer  was  assumed  to  occur  at 
both  ends  of  the  test  tube  (the  distance  Lu ) .  This  part  of 
the  tube  was  considered  to  be  a  straight  fin  of  uniform 
cross  section.  For  this  cylindrical  geometry: 


A   =  tt  (D*  -  D2)/4 


(4) 
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A,  =   P  Lu  (5) 


s     r   u 


where : 


A   =   cross- sect ional  area  of  tube  (m2 ) 

c  v  ' 

D   =   tube  outside  diameter  (m2 ) 

o  v  ' 

D,  =   tube  inside  diameter  (m2 ) 

A   =   tube  outside  surface  area  of 

s 

nonboiling  section  (m2 ) 
p   =   tube  outside  wall  perimeter  (m) 

P   =  -n-  D0 

Lu  =   non-boiling  length  of  the  test  tube  (m) 

It  was  assumed  that  the  temperature  at  the  base  of  the 
straight  fin  was  equal  to  the  average  wall  temperature  of 
the  active  boiling  section.  The  average  outer  wall  tempera- 
ture of  the  active  boiling  section  was  calculated  using  the 
radial  conduction  equation  from  the  thermocouple  position  to 
the  surface  of  the  tube. 

8 

Tavg=    (  X>)/8  (6) 

n=l 

Two  =   Tavg-  (^(  ln(l^/  Dx)/(2  ^  L  kc)  )    (7) 


Tj-,     Two 


where : 

1^  =   temperature  of  the  thermocouple 
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location  (°C) 

Tav=   average  wall  temperature  at  the 
thermocouple  locations  (°C) 

D-,  =   diameter  at  the  position  of  the 
thermocouple  (m) 

D2  =   outer  diameter  of  the  boiling  tube  (m) 

L   =   active  boiling  tube  length  (m) 

kc  =   thermal  conduct ivitiy  of  the  copper  (W/m.K) 

iL0  =   outer  wall  temperature  of  the  boiling 
test  tube  (°C) 

Tb  =   temperature  at  the  base  of  straight  fin  (°C) 

The  temperature  distribution  of  a  fin  of  uniform  cross- 
sectional  area  is  given  by  [Ref.  20]  with  the  following 
assumptions : 

1.  One-dimensional  conduction  in  the  x  direction. 

2.  Steady-state  condition. 

3.  Constant  thermal  conductivity. 

4.  Negligible  radiation  from  the  surface  of  the  fin. 

5.  No  heat  generation  in  the  fin. 

6.  Uniform  convective  heat  -  transfer  coefficient  over  the 
fin  surface. 

d28  /  dx2   -   m2e  =   0  (8) 


where 


0(x)  -   T(x)   -   Tsat  (9) 
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and 


m2   =   h  p  /  kcAc  (10) 


where : 


T__j=   saturation  temperature  of  fluid  (°C) 

h   =   convective  heat  -  transfer 
coefficient  (W/m2.K) 


Now,   assuming  negligible   heat  loss  from  the   thin  tube 
tip  (at  x  =  Lu) 

(d9  /  dx^=L=   0  (11) 


Ju 


and  also  at  the  base 


I      =       -k  A  (dT/dx)   =   -koAc(d0/dx)  (12) 

a  c  c        x=0      °  c         x=0 


where 


QB  =   heat- transfer  rate  from  the  base  (W) 

Using  the  above  boundary  conditions  (Eq.  11,  12),  the 
temperature  distribution  along  the  straight  fin  can  be 
expressed  as : 

(9/9b)  =  (  Cosh(m(Lu-  x))  /  Cosh(mL  )  )         (13) 


53 


and  the  fin  heat  -  transfer  rate  is: 

Qf  =   Qb  =  -Jh   P  ^cAc   9bTanh(mLu)  (14) 

where : 

QF  =   heat  -  transfer  rate  through  one  non-boiling 
length  of  the  test  tube  (W) 

The  average  difference  between  the  wall  temperature  and 
the  liquid  temperature  may  now  be  determined  from  the 
following  equation: 

(0/0b)  =   (|(Cosh(m(Lu-  x))/Cosh(mLu)  )  dx  )  /  L  u    (15) 

9   =   (0b/  mLu)  Tanh(mLu) 

9      Two  -   Ts  a^- 

where: 

T^0  -   T3at=   average  difference  between  outer 

wall  temperature  and  liquid  saturation 
temperature  (K) 

A  free  convection  correlation  stated  by  Churchill  and 
Chu  [Ref.  20]  was  applied  for  the  average  Nusselt  number  on 
a  horizontal  cylinder: 


_    V6 

I  0.387    Ra  \    2 

Nun         =      -j  0.60+ Q°         o/n-7  (16) 

Do  |  [1  +(0.559/Pr)9/16]8/27( 


-5         —  12 

10  <    RaD      <    10 

o 
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The  average  Nusselt  number  is 


h  Do  (17) 

NUD   =  "IT" 
o 


where : 


k   =   thermal  conductivity  of  R-114  (W/m.K) 


Solving    for   h 


|  0.387    Raj/6  j  2  (18) 

h      =      -r—    >  0  .  60  +  - 


Do   I     '  U+(0.559/Pr)9/16]8/27) 


Pr      =      * 

a 


gB(Tu.-TeM1)D 
Ra 


W0       SATy     o 


D  v    a 

o 


ft  1    Ap 

6     ~     "  p  It 


where 


Pr  =   Prandtl  number 

v   =   kinematic  viscosity  (m2/s) 

a   =   thermal  diffusivity  (m2/s) 


55 


(19) 


(20) 


(21) 


RaDo  =   average  Rayleigh  number 

gravitational  acceleration  (m/s2) 

volumetric  thermal  expansion 
coefficient  (1/K). 

density  of  R-114  liquid  (kg/m3) 


Do 

g 
3 


Now,    substitution   of    Equations    10,     15,       20    into    Equation 
18    results    in: 


-.3 „t   hP./Z^pi/e 


qSD    P    TanhAr- 

V         O     B  (   K 


c"c 


h      =      ~  <  0.60  +  0.3E 

o 


vaLu(^-s-) 


[l+(0.559/Pr)9/16]8/27 


(22) 


Equation  22  is  solved  for  h  by  an  iterative  technique 
within  a  range  of  precision  of  0.001.  Knowing  the  value  of 
the  natural  convection  heat -transfer  coefficient  along  the 
non-boiling  ends,  the  total  heat-loss  rate  is  calculated 
from  Equation  14,  and  the  heat  -  transfer  rate  through  the 
boiling  section  is  obtained  by  subtracting  the  total  heat- 
loss  rate  from  the  total  heat  -  transfer  rate: 


\ 


OSS 


2x 


% 


(23) 


Q   -   V  2xQF 


(24) 


56 


Finally,  the  heat  flux  from  the  boiling  surface  is: 


q    Q  /  a. 


(25) 


where 


Ab  =   7TDQL 
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V.  RESULTS  AND  DISCUSSION 

A.  OUTLINE  OF  THE  DATA  RUNS 

Using  the  procedure  outlined  in  Chapter  3,  nine  data 
runs  were  completed,  primarily  to  de-bug  the  experimental 
apparatus  for  its  successful  operation.  For  all  data  runs, 
the  same  smooth  copper  tube,  described  in  Chapter  2,  was 
used  as  the  boiling  surface.  Each  data  run  consisted  of 
seven  different  heat  fluxes  (1000,  2000,  5000,  10000,  20000, 
50000  and  95000  W/m2 )  with  a  specified  saturation  tempera- 
ture (0,  10  or  20  °C).  All  the  raw  data  were  stored  in  a 
user- specif ied  file,  named  "WHxx" ,  where  "W"  indicates 
Wieland  tube,  "H"  indicates  hard  copper  and  "xx"  indicates 
the  data  run  number.  Also,  the  heat- transfer  coefficient 
versus  heat  flux  data  were  stored  in  a  plot  file  named 
"Pxx".   A  summary  of  these  data  runs  is  given  in  Table  II. 

B.  LONGITUDINAL  AND  CIRCUMFERENTIAL  TEMPERATURE  VARIATIONS 

During  all  data  runs,  especially  at  high  heat  fluxes 
(greater  than  10  kW/m2 ) ,  large  temperature  variations  were 
observed  along  the  active  boiling  length  of  the  test  tube. 
Figure  5.1  shows  the  measured  wall  temperatures  at  each  of 
the  thermocouple  locations,  while  Figure  5.2  shows  the 
temperature  variations  in  the  circumferential  direction  (at 
the  axial  mid  point)  at  two  different  heat  fluxes  (2  and  92 
kW/m2 ) .  Despite  the  uniform  heat  flux  provided  by  the 
heater,  the  temperature  measurements  showed  considerable 
variations  in  both  axial  and  circumferential  directions. 
The  reason  for  these  temperature  variations  may  be  explained 
as  follows:  As  explained  in  Chapter  2,  during  the  manufac- 
turing  of  the   boiling  test   tube,   the   copper  sleeve   was 
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inserted  into  the  boiling  test  tube  without  any  shrink- fit 
process.  The  slip-fit  process  (with  a  diameter  clearance  of 
about  0.01  mm)  used  in  this  study  was  believed  to  have 
caused  considerable  thermal  contact  resistance  at  the  inter- 
face between  the  boiling  tube  and  the  copper  sleeve.  Also, 
it  is  quite  possible  that  this  contact  resistance  was  not 
uniform  along  both  the  axial  and  circumferential  directions. 
The  circumferential  temperature  distribution  may  also  be 
attributed  to  a  circumferential  variation  in  boiling  heat- 
transfer  coefficient.  During  the  high  heat  flux  boiling 
regime,  it  was  observed  that  the  rate  of  bubble  formation  at 
the  top  of  the  horizontal  boiling  tube  was  lower  than  at  the 
rest  of  the  boiling  tube;  i.e.,  the  boiling  heat  -  transfer 
coefficient  at  top  of  the  boiling  tube  was  relatively  lower 
than  the  other  locations.  This  poorer  performance  at  the  top 
of  the  tube  may  be  attributed  to  the  fact  that  this  portion 
was  not  receiving  sufficient  amount  of  fresh  liquid  to 
replace  the  liquid  that  evaporated.  The  obstruction  of  the 
liquid  flow  could  be  caused  by  the  growing  vapor  blanket 
from  the  bottom  to  the  top  of  the  tube.  Based  on  this  obser- 
vation, it  is  expected  that  the  temperature  at  top  of  the 
boiling  tube  could  be  greater  than  the  other  locations.  As 
shown  in  Figure  5.2,  location  5,  which  represents  the  top  of 
the  horizontal  boiling  tube,  measured  the  highest  tempera- 
ture difference  as  expected  from  the  above-mentioned  obser- 
vation . 

However,  positions  1  and  2  showed  a  sizable  discrepancy 
at  high  heat  flux,  and  this  was  due  most  likely  to  the  non- 
uniform contact  resistance.  Also,  a  comparison  of  thermocou- 
ples .2,  7  and  8  show  a  similar  discrepancy.  It  is  clear 
therefore,  that  the  shrink-fit  process  must  be  used  in  the 
future . 


59 


C.   PLOT  ANALYSIS  OF  NUCLEATE  BOILING  REGIME 

Figure  5.3  shows  a  typical  nucleate  pool-boiling 
performance  curve  of  the  smooth  copper  tube  in  R-114.  The 
observed  behavior  of  this  process  is  analyzed  by  studying 
the  different  regions  of  the  boiling  curve. 

From  point  A  to  point  B,  a  continuous  increase  in  wall 
superheat  (\jQ  -^at)  is  observed  when  heat  flux  is  increased. 
Also,  during  the  experimental  runs,  no  bubbles  were  observed 
along  the  boiling  tube  in  this  region  of  the  curve.  This 
region  corresponds  to  the  natural- convection  process. 

From  point  B  to  point  C,  a  reduction  in  wall  superheat 
is  observed,  while  the  heat  flux  continuously  increased. 
While  the  middle  portion  of  boiling  tube  showed  bubbles,  the 
rest  of  the  boiling  tube  showed  no  bubbles  during  the  exper- 
iment in  this  region.  This  region  is  known  as  the  mixed 
boiling  region,  where  transition  from  natural  convection  to 
nucleate  pool  boiling  heat  transfer  takes  place.  The  wall 
superheat  continued  to  decrease  until  numerous  nucleation 
sites  became  active  (point  C). 

After  point  C,  the  wall  superheat  starts  to  increase 
with  increasing  heat  flux  as  shown  by  region  C  to  D.  A  very 
high  density  of  bubble  formation  was  observed  along  the 
active  boiling  section  of  the  test  tube  (only  a  few  bubbles 
were  observed  on  the  nonboiling  sections  of  the  test  tube). 
See  Figure  5.4.  This  observation  confirmed  the  validity  of 
the  natural  convection  model  assumed  for  data  reduction. 
However,  the  R-114  liquid  returning  to  the  evaporator  from 
the  condenser  caused  a  mild  turbulance  in  the  entrance 
region  of  this  return  line  in  the  evaporator.  Inherently, 
this  affects  the  natural- convective  heat  -  trans fer  perform- 
ance of  nonboiling  ends  of  the  test  tube. 

When  heat  flux  is  gradually  decreased,  the  curve  follows 
a  different  path  (from  D  to  E)  as  shown  in  Figure  5.3.   This 
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is  due  to  the  existing,  stable  nucleation  sites  remaining 
active  for  a  wide  range  of  heat  fluxes.  Similar  analysis  can 
also  explain  Figure  5.5,  which  represents  the  same  data  run 
but  on  a  different  basis  (heat  -  transfer  coefficient  versus 
heat  flux) . 

D.  REPRODUCIBILITY  TEST  OF  THE  APPARATUS 

In  order  to  test  the  reproducibility  of  the  experimental 
apparatus ,  two  runs  were  performed  on  two  different  days 
(data  run  number  5  and  7  in  Table  II)  at  the  same  saturation 
temperature.  Figure  5.6  and  Figure  5.7  show  the  comparison 
of  these  two  data  runs .  It  can  be  seen  from  these  figures , 
that  there  is  very  good  agreement  between  these  two 
different  runs.  This  agreement  shows  the  ability  of  the 
apparatus  to  reproduce  data  runs  revealing  successful 
operation. 

E.  BOILING  PERFORMANCE  OF  SMOOTH  COPPER  TUBE  IN  R-114 

In  order  to  test  the  validity  of  the  data  taken  from  the 
present  experimental  apparatus  (despite  the  presence  of 
errors  owing  to  the  thermal  contact  resistance  mentioned 
earlier) ,  an  attempt  was  made  to  compare  the  present  data 
with  data  found  in  the  literature.  For  this  purpose,  data 
run  number  5  (10  °C  saturation  temperature  with  decreasing 
heat-flux  rate)  was  compared  with  two  sources  in  the  litera- 
ture as  discussed  below. 

1.   Comparison  with  Chongrungreong- Sauer  Correlation 

A  correlation  proposed  by  Chongrungreong  and  Sauer 
[Ref.  7]  for  the  nucleate  boiling  performance  of  refrige- 
rants and  ref rigerant- oil  mixtures  is  compared  with  run 
number  5.  This  correlation  is  based  on  a  dimensional  anal- 
ysis and  data  from  various  sources:  [Ref.  7],  [Ref.  14], 
[Ref.  15],  [Ref.  22]  and  [Ref.  23]. 
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where 


h  =  heat- transfer  coefficient  (W/m2.K) 
Q  =  heat  rate  (W) 
A  =  surface  area  (in2  ) 

D  =  characteristic  dimension  of  heated  surface  (m2 ) 
yT  =  viscosity  of  saturated  liquid  (gr/m.s) 

Li 

fg  =  latent  heat  of  vaporization  (W.s/g) 
C  =  specific  heat  of  liquid  (J/g.K) 
kL=  thermal  conductivity  of  liquid  (W/m.K) 
P  =  pressure  (atmospheres) 


=  volume  fraction  of  pure  refrigerant 


Pl  =  liquid  density  (g/cm3) 


P   =  vapor  density  (g/cm3) 


Eq ' n  (26)  agrees  very  well  with  theirdata  [Ref.  7] 
for  R-ll  as  well  as  with  Stephan's  data  [Ref.  22]  for  oil- 
free  refrigerants,  R-ll,  R- 12 ,  R-13,  R-21,  R-22,  R-113  and 
R-114  with  errors  smaller  than  16  percent.  Figure  5.8  pres- 
ents a  comparison  of  the  current  experimental  data  with  eq ' n 
(26).  As  can  be  seen  from  this  figure,  the  current  data  are 
in  excellent  agreement  with  this  predictive  equation.  (The 
analysis  given  in  Appendix  I  shows  that  the  uncertainty  in 
the  boiling  heat- transfer  coefficient  for  this  experimental 
investigation  is   about  20%.)    This  agreement   was  somewhat 
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unexpected  since  the  present  data  includes  some  kind  of 
contact  resistance  in  it.  If  this  contact  resistance  were 
known,  and  were  removed  from  the  data,  a  higher  boiling 
coefficient  would  result. 

2 .   Comparison  with  Data  of  Henrici 

The  nucleate  boiling  heat- transfer  coefficient  of 
R-114  and  R-114-oil  mixtures  from  a  horizontal,  smooth 
copper  tube  (30  mm  in  outer  diameter)  has  been  experimen- 
tally investigated  by  Henrici  [Ref.  10].  Since  Henrici's 
data  were  taken  from  a  tube  with  a  larger  diameter  than  the 
current  boiling  tube  diameter,  his  data  were  revised  to 
include  this  difference.  For  this  purpose,  a  semi- emprical 
correlation  to  show  the  effect  of  tube  diameter,  as  devel- 
oped by  Cornwell  et  al.   [Ref.  21],  was  used: 

Nu  =  CRe2'3  (27) 

where : 

Re  =  (q  D)  /  (hf gV )     and    Nu  =  (h  D)  /  k 
C  =  150  for  refrigerant 

Using  this  relationship,  Henrici's  data  were  revised 
for  a  diameter  of  15.875  mm  (5/8  in.)  used  in  this 
experiment . 

h(D=15.875  mm)  =  h(D=30  mm)  ( 30/ 15 . 875 ) l  ' 3 
h(D=15.875  mm)  =  1.24  h(D=30  mm) 
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A  line  representing  the  revised  Henrici's  data  is 
plotted  in  Figure  5.9  together  with  the  unrevised  data  line. 
It  can  be  seen  that  the  boiling  heat  -  transfer  coefficient 
increases  by  about  25  percent  when  the  diameter  is  decreased 
from  30  mm  to  15.875  mm  (50-percent  reduction  in  diameter). 

Figure  5.10  provides  Henrici's  natural  convection 
data  and  Henrici's  revised  boiling  data  in  comparison  with 
the  present  experimental  data.  As  expected,  the  natural- 
convective  heat- transfer  coefficient  is  lower  than  the 
boiling  heat- transfer  coefficient.  The  boiling  heat- transfer 
coefficient  from  Henrici's  data  is  higher  than  the  coeffi- 
cient found  in  this  study.  The  reason  for  this  may  be 
explained  as  follows.  In  Henrici's  experimental  investiga- 
tion, the  measurement  of  the  surface  temperature  has  been 
done  using  thermocouples  outside  of  the  boiling  surface. 
Unfortunately,  he  did  not  describe  in  detail  how  he  measured 
the  wall  temperatures.  If  the  thermocouples  were  not 
totally  imbedded  in  the  wall,  they  may  produce  incorrect 
wall  temperature  values.  If  a  portion  of  a  thermocouple  is 
in  direct  contact  with  the  liquid,  this  thermocouple  will 
measure  a  value  lower  than  the  actual  wall  temperature  and 
such  an  error  would  decrease  the  wall  superheat,  thus 
resulting  in  a  higher  heat- transfer  coefficient.  Also,  if 
the  contact  resistance  effect  is  removed  from  the  present 
data,  then  the  heat  transfer  coefficient  would  increase 
toward  the  Henrici  data. 

F.   EFFECT  OF  PRESSURE 

The  effect  of  pressure  on  the  nucleate-boiling  heat- 
transfer  performance  of  R-114  was  investigated  for  reduced 
pressures  (P/Pc  )  of  0.0268,  0.0391  and  0.0556  (data  run 
numbers  5,  8  and  9,  respectively  in  Table  II).  Figure  5.11 
shows   the    relative   increase   in    boiling   heat- transfer 
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coefficient  with  increase  of  pressure  as  expected.  In  order 
to  calculate  the  pressure  effect  on  the  boiling  performance, 
a  simple  equation  developed  by  [Ref .  7]  was  used. 


where 


h  =  6.17  (q) 


a.  5  5 


(</>) 


3.6  5 


0.2  k 


(28) 


h  =  heat- transfer  coefficient  (W/m2.K) 
q  =  heat  flux  (W/m2 ) 

0  =  volume  fraction  of  pure  refrigerant 
P  =  boiling  pressure  (atmospheres) 
Using  this  equation,   the  pressure  effect  on  the  boiling 
heat  transfer  coefficient  can  be  written  as  follows: 


h(Px)  =  h(P2)  (P1/P2  ) 


0.2  It 


According  to  the  above  equation,  the  effect  of  pressure 
from  127  kN/m2  to  87  kN/m2  decreases  the  boiling  heat- 
transfer  coefficient  by  about  ten  percent,  which  agrees 
reasonably  well  with  the  current  experimental  data. 
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Figure  5.2   Circumferential  Temperature  Variation 

on  the  Boiling  Tube. 
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VI.  CONCLUSIONS 

An  experimental  apparatus  to  study  the  nucleate- 
boiling  heat- transfer  performance  of  R-114  from  a 
single  horizontal  tube  has  been  designed,  constructed 
and  instrumented. 

An  analysis  scheme  to  establish  the  nucleate  boiling 
heat- transfer  coefficient  of  R-114  has  been  outlined, 
and  data  acquisition/reduction  programs  have  been 
written  and  tested. 

The  experimental  apparatus  was  successfully  operated 
in  the  saturation  temperature  range  of  0  °C  to  20  °C 
and  in  the  heat  flux  range  of  1  to  95  kW/m2 . 
The  test  runs  performed  on  different  days  on  the 
apparatus  showed  repeatabilty  within  ±  1  percent, 
revealing  successful  operation  of  the  apparatus. 
The  uncertainty  in  wall  superhet  (about  20%)  drasti- 
cally affected  the  boiling  heat- transfer  coefficient. 
This  can  be  attributed  to  the  large  temperature  vari- 
ations measured  along  and  around  the  boiling  tube 
because  of  a  non-uniform  contact  pressure  between  the 
copper  sleeve  and  the  boiling  tube. 

The  data  taken  on  this  experimental  apparatus  have 
been  compared  with  the  experimental  correlation  for 
refrigerant  and  refrigerant-oil  mixtures  developed  by 
Chongrungreong  and  Sauer  [Ref.  7].  The  present  data 
agreed  to  within  5  percent  with  this  correlation. 
The  current  data  were  also  compared  with  Henrici's 
[Ref.  10]  experimental  nucleate  boiling  data  (from 
smooth  horizontal  copper  tube  in  R-114).  The  present 
data  lie  as  much  as  50  percent  below  Henrici's  data. 
This  disagreement  may  be   attributable  to  the  thermal 
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contact  resistance  present  in  the  boiling  tube,  or  to 
the  errors  in  Henrici's  wall  temperature 
measurements . 
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VII.  RECOMMENDATIONS 

To  prevent  considerable  temperature  variations  along 
and  around  the  boiling  tube,  a  shrink-fit  process 
must  be  applied  in  order  to  have  more  uniform  contact 
pressure  between  the  copper  sleeve  and  the  boiling 
tube  . 

To  determine  the  effect  of  contact  resistance 
(between  copper  sleeve  and  boiling  tube)  on  the 
nucleate-boiling  heat  -  transfer  performance  of  a  hori- 
zontal tube,  a  series  of  experimental  investigations 
should  be  performed  with  different  contact  pressures. 
To  maintain  constant  electrical  power  input  into  the 
boiling  tube,  a  voltage  regulator  should  be  added  to 
the  input  line. 

To  keep  the  temperature  of  the  water- ethylene  glycol 
mixture  at  a  constant  value,  the  1/2-ton  capacity 
R-12  refrigeration  unit  should  be  replaced  with  a 
2- ton  unit . 

The  positive-displacement  pump  should  be  replaced 
with  a  centrifugal- type  pump,  in  order  to  deliver 
lower  flow  rates  without  reaching  an  electrical  over- 
load limit . 

To  check  the  saturation  temperature  of  the  boiling 
liquid  corresponding  to  the  measured  saturation  pres- 
sure, a  more  accurate  pressure  transducer  should  be 
used . 

The  active  length  of  the  cartridge  heater  should  be 
carefully  measured  or  obtained  from  the  manufacturer. 
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APPENDIX   A 
PHYSICAL    PROPERTIES    OF    FREON    FLUOROCARBON    COMPOUNDS 


"FREON"  11 

"FREON"   12 

"FREON"  113 

"FREON"   114 

Chemical  Formula 

CCI.F 

CC!  F, 

CCIF  CCIF, 

CCIF,  CCIF, 

Molecular  Weight 

13737 

120  92 

187  38 

170  93 

Boiling  Point  it  1  it™                                 *C 

"F 

23  82 
74  87 

-29  79 
-21  62 

4757 
117  63 

3.77 
38  78 

Freertng  Point                                                 *C 

•F 

-111 

-  168 

-  158 
-252 

-35 
-31 

-94 
-137 

Critical  Temp*ratur»                                       "C 

•F 

198  0 
388  4 

1120 
233  6 

214  1 
4174 

145  7 
294  3 

Critical  Pressure                                         atm 

lb!    tq  In   I'll 

435 

639  5 

40  6 
596  9 

33  7 

495 

32  2 

4732 

Critical  volume                                     ec  /mot 
co  ft/lb 

247 
0  0289 

217 

0  0287 

325 

0  0278 

293 

0  02  75 

Critical  Density                                         |/ce 
Ibs/cu  tt 

0  554 
34  S 

0  558 
34  8 

0  576 
36  0 

0  582 
36  32 

Density.  Liquid                                          g/cc 

itzyc(irr)                    rbs'cu  tt 

at  54  5*C  (130'F)                                g/cc 
Ibs/cu  It 

1476 
92  14 
1  40 
87  40 

1311 
81  84 
1.19 

74  28 

1  565 
97  69 
1  49 

93  01 

1  456 
90  91 
1  36 

84  90 

Density.  Sat'd  Vapor                                   g/| 
at  Boiling  Pom)                           Ibs/cu  tt 

5  86 
0367 

6  33 
0  395 

7  38 
0461 

7  83 
0  489 

Specific  Heat,  liquid 
(Heat  Capacity)                       eal/(g)   (*C) 
at  25-C  (77'F)                      Btu/(lb)  ('Fl 

0  208 

0732 

0218 

0  243 

Specific  Heit,  Vapor, 

at  Const   Pressure   (1   atm)        <>'!()    CO 

at  25'C  (77'F)                      Btu/(lb)  CO 

0142  •    (iw-n 

0  145 

0161  "(Mo-n 

0  170 

Specilic  Heat  Ratio  ol  Vapor 

at  25'C  and  1  atm                           Cp/Cv 

1 137 »   (13^;cn 

1  137 

'■""•(Mo'-n 

1  084 

Heal  ot  Vaporisation                                eal/g 
at  Boiling  Point                                 Btu/lb 

43  10 
7751 

39  47 

71  04 

35  07 
63  12 

32  51 

58  53 

Thermal  Conductivity  at  25'C  (77'F) 
Btu/(hr)  (ft)('F) 
Liquid' 

Vapor  (1  aim) 
(Data  Irom  ASHRAE  In  most  cases) 

0  0506 
0.00451 

0  0405 
0.00576 

0  0434 

0C')14 

(0  5  atm) 

0  0372 
0  0060 

Viscosity  at  25'C  (77'F) 

Liquid                                           centipoisa 
Vapor  (I  atm)                              centipoisa 
(Data  liom  ASHRAE  In  most  cases) 

0  430 
00107 

0214 
0  0123 

0.68 

onio 

(0  1  aim) 

0  36 
00112 

Surlace  Tension  at  25'C  (77'F) 
dynes/cm 

IS 

9 

173 

12      » 

Refractive  inrVi  ol  Liquid 
at  25'C  (77'F) 

1  374 

1  287 

1  354 

1  288 

Relative  Dielectric  Strength  of  Vapor  at  1  atm 
and  25'C  (77'F)  (nitrogen  =  1) 

3  71 

246 

3  9  (0  44  atm) 

3  34 

Dielectric  Constant 
Liquid 
Vapor  (1  atm)"« 

2  28  9  29'C 
1  0036  <9  24*C* 

1^32          ®29'C 

2  41   9  25"C 

2  25  0  25"C 
1.0043  a  26  8"C 

SoluMllty  ol  "Freon"  in  Water 
at  1  atm  and  25'C  (77'F)                art  % 

0.11 

0  028 

0017 
(Sat'n  Pres) 

0013 

Soluhility  ol  Water  in  "FreoiT 
at  25'C  (77'F)                                   *t  % 

0  011 

0  009 

0011 

0  009 

Solubility  Parameter  (    ) 
Kauri  Bulanol  Value  (KB) 

75 
60 

61 
18 

72 

32 

62 
12 

Toaicity"                                           ppm  (y/i) 
Threshold  limit 

Value  (TLV)                                    m|/in' 

1000 
5600 

1000 
4950 

1000 
7600 

1000 
7000 
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APPENDIX   C 
APPLICATIONS    OF    "FREON"    FLUOROCARBON    COMPOUNDS 


The  table  t>elow  is  intended  to  provide  3  general  view  of  the  range  of  applications  and  is 
For  specialized  applications  or  more  detail,  please  make  specific  request. 

not  all  inclusive. 

Fluor  ocarbon 

Refrigerants 

Aerosol  Propellantslbl 

Solvents    Blowing  Agents. 

Fire  E  xtinguishantj. 

Dielectric  Fluids  and  Other  Uses 

"Freon"  14 







"Frenn"  23 

Component  of  "Fteon"  503 
a.'eoitope13' 

— 



"Freon"  13 

Specialty  low  temperature 
applications,  usually  in  "cascade" 
systems 

— 



"Freon"  116 



Dielectric  fluid 

"Freon"  I3BI 

Intermediate  between  "Freon"  13 
and  "Fteon"  22  lor  medium  to  low 
temperature  applications.  Not 
extensively  used 

— 

Efficient  fire  extmguishant 
(Halon'  13011  especially  suited  for 
aulomatic  protection  ol  materials 
subiect  to  water  damage  and  ot  areas 
occupied  by  personnel 

"Freon"  22 

Household  and  commercial 
tefriqetation  and  air  conditioning 
applications   Permits  use  ol  smaller 
equipment.  Component  of 
azeotropes,dl 





Freon"  1 15 

Used  as  an  azeofrope  component  in 
"Fieon  '  502la'. 

Accepted  as  a  lood  propellant  by  the 
FDA.  this  material  is  well  suited  for 
food  aerosols  and  linds  use  in  lat 
emulsion  food  whips   Good  loam 
stability  with  absence  ol  odor  or  taste 

Oielectric  fluid,  an  economic 
teplaceinent  for  "Fteon"  116  in 
most  dielectric  applications. 

"Freon"  12 

Most  widely  used  refrigerant  in 
household,  automotive  and 
commercial  refrigeration  and  air 
conditioning  systems  Also  as  a 
coniponent  ol  a/co^pes!*1  and, 
in  high  purity  form  ("Freon" 
Ireezant]  approved  as  a  direct 
contact  freezing  agent  for  foods. 

Most  widely  used  high  pressure 
propellant  lot  non  food  use   8lends 
with  "Freon"  11  and  "Fteon"  114 
are  widely  used 

8lowing  agent  lot  loamed  plastics 
applications!''!    Oielectric  gas 

"Freon"  1 14 

In  large  industrial  process  cooling 
and  air  conditioning  systems  using 
multi  stage  centrifugal  compressor*. 

Low  pressure  propellant  alternative 
to  "Freon"  1  1,  having  poorer 
solubility  properties  and  less  odor. 
Especially  used  in  personal  products 

Blowing  agent  lor  loamed  plastics. 

"Freon"  21 



Heat  transfer  fluid 

"Freon"  1 1 

Widely  used  in  centrifugal 
compressors  lor  industrial  and 
commercial  air  conditioning  systems 
and  for  industrial  cooling  of  process 
water  or  brine   Low  viscosity  and 
freezing  point  permit  use  as  a 
low  temperature  cooling  liquid. 

Most  widely  used  low  pressure 
propellant  for  non  lood  use   Does 
not  provide  adequate  piessure  alone, 
so  is  almost  entirely  used  in  blends 
with  "Freon"  12(c). 

Occasionally  used  as  a  solvent 
("Freon"  MF)    Blowing  agent  lor 
foamed  plastics' r" 

"Freon"  113 

In  commetcial  and  industrial  air 
conditioning  and  process  water  or 
brine  chilling  using  cenrrtlugal 
compressors,  par ticularfy  in  small 
tonnage  applications. 

Solvent  in  some  aerosol  lormulations, 
usually  propelled  with  "Freon"  12 

Extensively  as  a  solvent  ("Freon"  TF) 
alone  and  in  special  purpose 
lormulations  for  a  wide  range  of 
critical  cleaning  needs  In  cutting 
fluid  formulations,  Valclene 
dry  cleaning,  etc 

(a)  A  numltfi  ol  azeobopes  ("Freon"  500.      Freon  '  50?    etc  1  are  available  tor  refrigeration  use    Bulletins  describing  lhe  composition, 
properties  an<1  uses  ol  these  mixtures  are  available  on  request. 

Mil  Norma'iv  aerosol  propellents  are  blended  to  give  the  required  vapor  pressure  and  solubility  requirements 

Id    'Freon  '  1 1  S  is  a  stabilised  grade  lienuently  used  in  formulations  conlatning  alcohols  or  water 

Idl  "F  reon     12  and  "Freon  '  M  ate  widely  used  as  blowing  agents  lot  a  range  of  foamed  plasties  in  which  they  ptovide  excellent  cell  structure, 
and  the  trapped  lluorocart>on  in  c  1  used  cell  loams  signidcanily  improves  insulation  and  recovery  properties  ot  the  loam 
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APPENDIX  D 
THERMOCOUPLE  CALIBRATION 

A.   EQUIPMENT  USED 

The  equipment  used  in  thermocouple  calibration  is  shown 
in  Figure  D.l,  and  a  brief  description  of  each  component  is 
given  below: 

1 .  Thermocouple  Wire 

Type-T  ( copper- constantan)  Teflon- coated  wire  of 
0.254  mm  (0.01  in.)  in  diameter  was  used  for  all  thermocou- 
ples . 

2 .  Calibration  Bath 

A  Thermos  flask  was  used  as  the  calibration  bath.  In 
order  to  maintain  an  isothermal  temperature  distribution,  a 
motor-driven  mixer  was  used.  Observations  during  the  cali- 
bration procedure  showed  bath  temperature  fluctuations  to  be 
±  0.002  K. 

3 .  Thermocouple  Readout 

A  Hewlett  Packard  3497A  automatic  data  acquisition/ 
control  system  and  a  Hewlett  Packard  9826  computer  were  used 
to  read,  analyze  and  record  the  calibration  data. 

4 .  Reference  Temperature 

A  Hewlett  Packard  2804A  quartz  thermometer  was  used 
to  measure  the  bath  temperature.  This  quartz  thermometer  had 
a  resolution  of  0.0001  K,  while  the  manufacturer-guaranteed 
accuracy  was  better  than  ±  0.03  K. 
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B.   PREPARATION  FOR  CALIBRATION 

1 .  Thermocouple  Preparation 

A  total  of  five  thermocouples  were  prepared  for  the 
calibration  run.  Two  thermocouples  were  made  from  the 
beginning  of  a  spool  of  copper- constantan  wire;  two  were 
made  from  the  end;  and  one  was  made  at  the  mid  section  using 
the  following  procedure: 

•  The  Teflon  insulation  was  removed  for  a  length  of  about 
4  mm  from  one  end  of  aim  long  piece  of  wire  and  a 
thermocouple  bead  was  made  using  a  Dynatech  Corporation 
thermocouple  welder. 

•  The  other  end  of  the  thermocouple  wire  was  connected  to 
the  data  acquisition  and  control  unit  through  a  junc- 
tion box. 

2 .  Computer  Program 

A  short  computer  program  (TCAL)  was  written  to 
accept  the  thermocouple  readings  through  the  data  acquisi- 
tion system  and  the  bath  temperature  through  the  digital 
quartz  thermometer.  A  listing  of  the  computer  program,  TCAL 
is  provided  In  Appendix  G.  This  program  prints  all  data  and 
the  discrepancy  (i.e.,  the  Quartz  thermometer  reading  minus 
the  thermocouple  reading)  as  well  as  it  stores  the  data  on  a 
computer  disk.  In  order  to  convert  the  e.m.f.  values  to 
temperature,  the  manufacturer's  conversion  equation  (D.l) 
was  used  in  the  computer  program  TCAL: 


T  =  a0+  a-^1  +  a2E2  +  a3E3  (D.l) 

+  a4E"  +  a5E5  +  a6E6+  a7E7 
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where 


T  =  temperature  (°C) 

ao=  0.100860910 

ax=  25727.94369 

a2  =  -767345.8295 

a3=  78025595.81 

a4=  -9247486589 

a5=  6.97688E+11 

a6  =  -2.66192E+13 

a7=  3.94078E+14 

E  =  thermocouple  reading  (volts) 


CALIBRATION  PROCEDURE 

•  Since  the  temperature  measurements  during  this  investi- 
gation ranged  from  25  °C  to  70  °C,  thermocouple  cali- 
bration was  performed  in  this  region.  First,  hot  water 
(at  about  75  °C)   was  added   to  the  thermos  flask  cali- 

•  bration  bath.  Following  proper  mixing,  all  thermocouple 
readings  were  recorded  through  the  data  acquisition 
system  and  the  bath  temperature  was  measured  by  the 
quartz  thermometer.  To  obtain  lower  temperatures, 
small  quantities  of  cold  water  were  gradually  added  to 
the  bath. 

•  The  discrepancy  (i.e.,  the  Quartz  thermometer  reading 
minus  the  thermocouple  reading)  was  plotted  against  the 
thermocouple  reading  as  shown  in  Figure  D.2.  The 
second-order  polynomial  curve  (D.2)  shown  in  this 
figure  was  generated  using  the  data  from  all  5 
thermocouples . 
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DCP  =  -6.7422934E-02  +  9 . 0277043E-03  T  -       (D.2) 
9.3259917E-05  T2 

where : 

T  =  thermocouple  reading 

(from  equation  D.l)   (°C) 
DCP  =  discrepancy  (K) 

Thus,   the   corrected  temperature   values  were   obtained 
using  equation  (D.3): 


DCP  +  T  (D.3) 


where 


T  =  actual  calculated  temperature  (°C) 
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APPENDIX  E 
DATA  REDUCTION  PROGRAM 


1090!  FILE  NAME:  DRPR1 

1010!  DATE:      October  19.  1984 

1020!  REVISED:    November  I.  1984 
1030! 

1040  BEEP 

1050  PRINTER  IS  1 

1060  PRINT  USING  "4X . ""Se 1 ec t  option: 

1070  PRINT  USING  "6X,""0  Taking  data  or  re-process inq  previous  data 

I H80  PRINT  USING  "6X,""1  Plotting  data 

1090  INPUT  Idp 

I  100  IF  ldp  =  0  THEN  CALL  Main 

1110  IF  Idp-1  THEN  CALL  Plot 

1120  END 

1  1  30  SUB  Main   - 

1  140  CUM  /Cc/  C(7) 

1 150  DIM  Emf < 12) .T( 12) ,D!a<2) .D2a<2) . Dia<2) ,Doa(2) . La<2) , Lua(2 ) ,Kcua<2> 

1160  DATA  0.10086091.25727.94369,-767345.8295,78025595.81 

1170  DATA  -9247486589.  6.  97688E+1  1  ,-?..  66I92E+13.  3.  94078E+14 

1  180  READ  C<«> 

! 190  PRINTER  IS  701 

1200  CLEAR  709 

1210  BEEP 

1220  INPUT  "ENTER  MONTH.  DATE  AND  TIME  (MM : DD :HH : MM :SS) " , DateS 

1230  OUTPUT  709:"TD":DateS 

1240  OUTPUT  70q:"TD" 

1250  ENTER  709;DateS 

1250  PRINT 

1270  PRINT 

1280  PRINT 

1290  PRINT  "  Month,  date  and  time  :":Date$ 

1300  PRINT 

1310  PRINT  USING  "1  OX  ,  ""NOTE :  Program  name  :  DRP 

1320  BEEP 

1330  INPUT  "ENTER  DISK  NUMBER". Dn 

1340  PRINT  USING  " 1 6X  ,  ""D isk  number   -  "",ZZ":Dn 

1350  BEEP 

1360  INPUT  "ENTER  INPUT  MODE  (0-3054A  .  1 -FILE)" . Im 

1370  IF  Im-0  THEN 

1  380  BEEP 

1330  INPUT  "GIVE  fl  NAME  FUR  THE  RAW  DATA  FILE" . D2_fi leS 

1400  PRINT  USING  "1 6X . ""New  file  name:  "" , 1 4A" : D2_f 1 leS 

1410  CREATE  BOAT  D2_fiIeS,20 

1420  ASSIGN  »File2  TO  D2_file$ 

I '130! 

1440!  DUMMY  FILE  UNTIL  Nrun  KNOWN 

14S0  D1_f i le$="DHMMY" 

1460  CREA1E  BDAT  D 1 _f i leS ,20 

1470  ASSIGN  @F  i le1  TO  D 1 _f  i leS 

1480  OUTPUT  «Fl lei :Date$ 

1490  BEEP 

1500  INPUT  "GIVE  A  NAME  FOR  THE  PLOT  FILE" . P_f i 1 e$ 

1510  CREATE  BDAT  P_f ileS,5 

1520  ASSIGN  SPIot  TO  P_fileS 

1530  BEEP 

1540  INPUT  "ENTER  NUMBER  OF  DEFECTIVE  TCS  < 0-DEFAULT )" . Id tc 

1550  IF  ldtc=0  THEN 

1560  Ldtc!=0 

1570  Ldtc2=0 

1580  PRINT  USING  "16X,""No  defective  TCs  exist 

1590  END  IF 
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1600  IF  Idtc=t  THEN 

1610  BEEP 

1620  INPUT  "ENTER  DEFECTIVE  TC  LOCATION" .Ldtc I 

1630  PRINT  USING  "16X,""TC  is  defective  at  location  "",D";Ldtc1 

1 640  Ldtc2=0 

1650  END  IF 

1660  IF  Idtc-2  THEN 

1670  BEEP 

1680  INPUT  "ENTER  DEFECTIVE  TC  LOCATIONS" . Ldtc 1 , Ldtc2 

1690  PRINT  USING  "16X.""TC  are  defective  at  locations  ""  .  D  .  4X  ,  D"  :  Ldtc 1  .  Ldtc2 

1700  END  IF 

1 710  IF  Idtc>2  THEN 

1720  BEEP 

1730  PRINTER    IS    1 

1740  BEEP 

1750  PRINT  "INVALID  ENTRY" 

1760  PRINTER  IS  701 

1770  GOTO  1530 

1 780  END  IF 

1790  OUTPUT  -3F  1 1  e  1  ;Ldtc1  .LdtcZ 

1800!  Im=1  option 

1810  ELSE 

1820  BEEP 

1830  INPUT  "GIVE  THE  NAME  OF  THE  EXISTING  DATA  FILE" ,D2_f i le$ 

1840  PRINT  USING  "16X.""IJld  file  name:  ""  .  1  4A"  :  D2_f  i  leS 

1850  ASSIGN  9Flle2  TO  D2_file$ 

1860  ENTER  8File2:Nrun 

1870  ENTER  W  i  le2 ; Dold$ 

1880  BEEP 

1890  INPUT  "GIVE  A  NAME  FOR  PLOT  F ILE" , P_f i leS 

1900  CREATE  BOAT  P_f ile$.5 

1 910  ASSIGN  ®Plot  TO  P_file'S 

!920  PRINT  USING  " 1 6X . ""This  data  set  taken  on  :  ""  ,  1  4A" : Do IdS 

1930  ENTER  «File2:Ldtc1  .Ldtc2 

I940  PRINT  USING  " 1 6X . ""Thermocoup les  were  defective  at  locat ions : "" . 2( 3D . 4X ) ' 
Ldtcl ,Ldtc2 

1950  ENTER  9File2:Itt 

1960  END  IF 

1970  PRINTER  IS  I 

1980  IF  lm=0  THEM 

'990  BEEP 

2000  PRINT  USING  "4X. ""Select  tube  type 

2010  PRINT  USING  "6X , ""0=Smooth  4  inch  Ref 

2020  PRINT  USING  "6X . "" 1 =Smoo th  4  inch  soft  Cu 

2030  PRINT  USING  "6X , ""2=Smoo th  8  inch  soft  Cu 

2040  INPUT  Itt 

2050  IF  Itt>2  THEN 

2060  BEEP 

2070  PRINT  "INVALID  ENTRY" 

2080  GOTO  2000 

2090  END  IF 

2100  OUTPUT  ©File!  :  Itt 

2110  END  IF 

2120  PRINTER  IS  701 

2130  PRINT  USING  "16X.""Tube  Type  is  :  "".D":Itt 

2140  BEEP 

2  I  50  INPUT  "ENTER  OUTPUT  VERSION  ( 0  =  LONG . 1 -SHORT >", lov 

2160! 

2170!  Dl-Diameter  at  thermocouple  positions 

2180  DATA  .01 1 125, .01 1 1 125, .01 143 

2190  READ  Dla(') 
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2200 

2210! 

2220! 

2230 

2240 

2250 

22S0! 

2270! 

2280 

2290 

2300 

2310! 

2320! 

2330 

2340 

2350 

2360! 

2370! 

2380 

2  330 

2400 

2410! 

2420! 

2430 

2440 

2450 

2460! 

2470! 

2430 

2490 

2500 

2510 

2520 

2530 

2540 


D1=D1a(Itt) 

D2=Diameter  of  test  section  to  the  base  of  fins 
DATA  .015875, .015875. .015875 
READ  D2a<*> 
D2=D2a(Itt) 

Di=Inside  diameter  of  unenhanced  ends 
DATA  .0127. .0127, .0127 
READ  Dia<*> 
Di=Dia<Itt> 

Do=0utside  diameter  of  unenhanced  ends 
DATA  .015875. .015875, .015875 
READ  Doa(*> 
Do=Doa(Itt) 


L=Length  of 
DATA  .1016, 
READ  La<*> 
L=La(Itt) 


enhanced  surface 
1016,  .2032 


Lu=Length  of  unenhanced 
DATA  .0254, .0254. .0762 
READ  Lua<«) 
Lu=Lua<Itt> 


surface  at  the  ends 


2580 

2590 

2600 

2610 

2620 

2630 

2640 

2650 

2660 

2670 

2680 

2690 

2700! 

2710! 

2720 

2730 

2740 

2  750 

2760 

2770 

2780 

2730 

2800 


Kcu=Thermal  Conductivity  of  tube 
DATA  398.344,344 
READ  Kcua(-) 
Kcu=Kcua(Itt> 
A=PI»(Do  2-Di  2)/4 
P=PI*Do 
J=1 
Sx  =  0 
Sy  =  0 
Sxs  =  0 
Sxy  =  0 
Repeat :  ! 
IF  lm=0  THEN 

ON  KEY  0.15  RECOVER  2530 
PRINTER  IS  1 
PRINT  USING  "4X 

USING  "9X 

USING  "9X 

USING  "9X 

USING  "4X 


PRINT 

PRINT 

PRINT 

PRINT 

BEEP 

INPUT 

IF  Idr 

LOOP 
IF  Id 
OUTPII 
BEEP 
INPUT 
PRINT 
Err  =  1 
FOR  I 
OUTPIJ 
Sum  =  0 


"SELECT  OPTION 

"0=TAKE  DATA 

"1=SET  HEAT  FLUX'" 

"2=SET  Tsat 

"NOTE 


KEY  0  =  ESCAPE 


Ido 

=0  THEN 


3630 


TO  SET  HEAT  FLUX 

o=1  THEN 

T  709:"AR  AF62  AL63  VR5" 

"ENTER  DESIRED  0dp".DqdP 

USING  "4X. ""DESIRED  Odp   ACTUAL  Odp" 
000 

=  1  TO  2 
T  709;"AS  SA" 
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2b!0 

FOR  Ji=1  TO  5 

2820 

ENTER  703:E 

2830 

Sum  =  Su.m  +  E 

2840 

NEXT  Ji 

2350 

IF  1=1  THEN  Volt-Sum»5 

2860 

IF  1=2  THEN  AmP=E 

2870 

NEXT  I 

2880 

fiqdp=Volt*Amp/(PI*D2*L) 

2890 

IF  ABS<Aqdp-DqdP>>Err  THEN 

2900 

IF  Aqdp>Dqdp  THEN 

2910 

BEEP  4000,  .2 

2320 

BEEP  4000.  .2 

2930 

BEEP  4000,  .2 

2940 

ELSE 

2950 

BEEP  250.  .2 

2960 

BEEP  250.  .2 

2970 

BEEP  250.  .2 

2980 

END  IF 

2990 

PRINT  USING  "4X.MZ.3DE.2X.MZ.3DE" 

3000 

WAIT  2 

3010 

GOTO  2780 

3020 

ELSE 

3030 

BEEP 

3040 

PRINT  USING  "4X.MZ.3DE.2X.MZ.3DE" 

3050 

Err=500 

3060 

WAIT  2 

3070 

GOTO  2780 

3080 

END  IF 

3090 

END  IF 

3100! 

3'  10! 

LOOP  TO  SET  Tsat 

3120 

IF  Ido=2  THEN 

3130 

OUTPUT  703;"AR  AF33  AL33  VR5" 

3140 

OUTPUT  709:  "AS  SA" 

3150 

BEEP 

3160 

INPUT  "ENTER  DESIRED  Tsat". Dt Id 

3170 

PRINT  USING  "4X.""D  Tsat    A  Tsat 

3180 

Old  =  0 

3190 

ENTER  709:E1 iq 

3200 

Atld  =  FNTvsv(El  iq) 

3210 

Psat=FNPsat<Atld> 

3220 

IF  ABS(Atld-Dtld)>.2  THEN 

3230 

IF  Atld>Dtld  THEN 

3240 

BEEP  4000.  ,2 

3250 

BEEP  4000.  .2 

3260 

BEEP  4000.  .2 

32  70 

ELSE 

3280 

BEEP  250.  .2 

3290 

BEEP  250.  .2 

3300 

BEEP  250.  .2 

3310 

END  IF 

3320 

Err=Atld-01d 

3330 

Old-Atld 

3340 

PRINT  USING  "4X.4(MDD.DD.3X)-";Dtl 

3350 

WAIT  2 

3360 

GOTO  3190 

3370 

ELSE 

3380 

IF  ABS<AtJd-Dtld)>. 1  THEN 

3330 

IF  Atld>Dtld  THEN 

3400 

BEEP  3000. .2 

3410 

BEEP  3000. .2 

:  Dqdp , Aqdp 


Dqdp  .Aqdp 


Diff     Psat" 


Atld.Err .Psat 
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3420 

ELSE 

3430 

BEEP  800.  .2 

3440 

BEEP  300,  .2 

3450 

END  IF 

34G0 

Err=Atld-Qld 

34  70 

Dld-fltld 

3480 

PRINT  USING  "4X.4<MDD.DD,3X)";Dtlc 

3490 

WAIT  2 

3500 

GOTO  3190 

3510 

ELSE 

3520 

BEEP 

3530 

Err=At Id-Old 

3540 

Old-Atld 

3550 

PRINT  USING  "4X,4<MDD.DD.3X>";Dtk 

3560 

WAIT  2 

3570 

GOTO  3190 

3580 

END  IP 

3590 

END  IF 

3600 

END  IF 

3610! 

ERROR  TRAP  FOR 

Ido  OUT  OF  BOUNDS 

3620 

IF  Ido>2  THEN 

3630 

BEEP 

3640 

GOTO  2620 

3650 

END  IF 

3660! 

3670! 

TAKE  DATA  IF  Im  = 

3680 

BEEP 

3690 

INPUT  "ENTER  BULK  OIL  7."  .  Bop 

3700 

OUTPUT  709:MAR 

AF25  AL36  VR5" 

3710 

FOR  1=1  TO  12 

3720 

OUTPUT  709: "AS 

SA" 

3730 

Sum=0 

3740 

FOR  Ji-1  TO  20 

3750 

ENTER  709:E 

3760 

Sum=Sum+E 

3770 

NEXT  Ji 

3780 

Emf <I)=Sum/20 

3790 

NEXT  I 

3800 

OUTPUT  709;"AR 

AF62  AL63  VR5" 

3810 

FOR  1=1  TO  2 

3820 

OUTPUT  709; "AS 

SA" 

3830 

Sum  =  0 

3840 

FOR  J 1-1  TO  20 

3850 

ENTER  709;E 

3860 

S  um  =  S  urn  +  E 

3870 

NEXT  Ji 

3880 

IF  1=1  THEN  Vr= 

Sum/20 

3890 

IF  1=2  THEN  Ir= 

■Sum/20 

3900 

NEXT  I 

3910 

ELSE 

3920 

ENTER  9File2:Bop.Told$.Emf <«> .Vr. I 

3930 

END  IF 

3940! 

3950! 

CONVERT  emf'S  TO 

3960 

Tua=0 

397  0 

FOR  1=1  TO  12 

3980 

IF  ldtc>0  THEM 

3990 

IF  I=Ldtc1  OR  I 

=Ldtc2  THEN 

4000 

T(I)=-99.99 

4010 

GOTO  4060 

4020 

END  IF 

Atld.Err ,Psat 


,Err ,Psat 
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■4UiU  LNL)  ih 

4040  T(I)=FNTvsv(Emf  (D) 

4050  IF  I<9  THEN  Twa= Tua+T< I ) 

4060  NEXT  I 

4070  Tua=Tuia/(8-Idtc) 

4080  Tld=T<9> 

4090  Tv=(T(10)+T(11)>/2 

4100  Tsump=T(12) 

4110  Amp=Ir 

4120  Volt=Vr*25 

4130  0=Volt«flmp 

4140  IF  Itt-0  THEN 

4150  Kcu.=  FNKcu(Tua) 

4  160  ELSE 

4170  Kcu=Kcua<Itt> 

4180  END  IF 

4190! 

4200!  FOURIER  CONDUCTION  EOUATION  WITH  CONTACT  RESISTANCE  NEGLECTED 

4210  Tw=Twa-Q-LQG<D2/D1 >/<2*PI*Kcu*L) 

4220  Thetab-Tu-Tld 

4230  IF  Thetab<0  THEN 

4240  BEEP 

4250  INPUT  "TWALL<TSAT  (0=CONTINUE.  1=END>".Iev 

4260  IF  lev=0  THEN  GOTO  2590 

4270  IF  Iev-1  THEN  5030 

4280  END  IF 

4290! 

4300!  COMPUTE  VARIOUS  PROPERTIES 

4310  Tf ilm-FNTf ilm(Tu,Tld) 

4320  Rho-FNRho(Tfilm) 

4330  Mu=FNMu(Tf ilm) 

4340  K-FNK(Tfilm) 

4350  CP=FNCp(Tf ilm) 

4360  Beta=FNBeta(Tf ilm) 

4370  Ni=Mu/Rho 

4300  Alpha=K/(Rho*Cp) 

4390  Pr=Ni/AlPha 

4400  Psat=FNPsat(Tld) 

4410! 

4420!  COMPUTE  HEAT  TRANSFER  COEFFICIENT 

4430  Hbar=i90 

4440  Fe=<Hbar~P/(Kcu»A>> ' .5*Lu 

4450  Tanh=FNTanh<Fe> 

4460  Theta=Thetab*Tanh/Fe 

4470!  PRINTER  IS  701 

4480!  PRINT  USING  "'4X . 7 ( 1 X , MZ . 3DE > " ; Hbar . Fe , Tanh . The  tab . The ta. Be ta, N i 

4490  Xx=(9.81»Beta»Thetab*Do" 3»Tanh/(Fe*Ni»Alpha>)  .166667 

4500  Yy=(1  +  < .559/Pr)  '  (9/ 1 6 > >  < 8/27 ) 

4510  Hbarc=K/Do«< . 6+ . 387-Xx/ Yy )  2 

4520  IF  ABS((Hbar-Hbarc)/Hbarc)>.001  THEN 

4530  Hbar=(Hbar+Hbarc)*.5 

4540  GOTO  4440 

4550  END  IF 

4560  Ql=(Hbar*P-Kcu*A)  . 5*The tab'Tanh 

4570  0c=0-2»01 

4580  As=PI-D2»L 

4590  QdP=Oc/As 

4600  Htabe=Qdp/Thetab 

4610! 

4620!  RECORD  TIME  OF  DATA  TAKING 

4630  IF  Im-0  THEN 
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4640  OUTPUT  7fl9;"TD" 

4650  ENTER  709;Told£ 

4660  END  IF 

4670! 

4680!  OUTPUT  DATA  TO  PRINTER 

4690  PRINTER  IS  701 

4700  PRINT 

4710  IF    lov=0    THEM 

4720  PRINT    USING    M  1  1  X  ,  ""Data   Set    Number    =    ""  ,  DDD  ,  2X  ,  ""Bulk    Oil    '/.    =    "",DD.D.5X.t 

4A": J.Bop.ToldS 

4730  PRINT 

4740  PRINT    USING    "MX.  ""TC   No:  1  2  3  4567 

8 

4750  PRINT  USING  " 1 ! X . ""Temp  : "" . 8< 1 X , MDD . DD > "  :  T <  I ) . T < 2 > . T<3 >  .  T < 4 >  ,  T < 5  )  .  T < 6  )  ,  T < 
7)  .T(8) 

4760  PRINT  USING  "MX,""  Tua    Tliqd    Thetab     Tvapr   Psat    Tsuitid 

4770  PRINT  USING  " 1 1 X , 2 (MDD . DD , 1 X ) . MZ . DDDE . 1 X , 2 ( 1 X , MDD . DD ) . 2X . MDD . D" ; Twa. T Id . Th 
e  tab , Tv . Psat , Tsump 

4780  PRINT  USING  "1IX,""  Htube      Qdp 

4790  PRINT  USING  "1 1 X ,MZ .3DE , 1 X . MZ . 3DE" : Htube . Qdp 

4800  END  IF 

4810  IF  Iov-1  THEM 

4820  IF  J- I  THEN 

4830  PRINT  USING  "MX.  ""RUN  No  Oil'/.   Tsat    Htube      Qdp        Thetab 

4840  END  IF 

4850  PRINT  USING  " 1 2X . DDD . 4X , DD , 2X . MDD . DD . 3d  X , MZ . DDDE) " ; J , Bop . T Id . Htube . Qdp . Th 

etab 

4  860  END  IF 

4870  IF  Im-0  THEM 

4880  BEEP 

4890  INPUT  "OK  TO  STORE  THIS  DATA  SET  ( !  =  Y . 0  =  N) ?"  .  Ok 

4900  END  IF 

4910  IF  Ok  - 1  OR  Im=l  THEN  J-J+l 

4920  IF  Ok  =  1  AND  lm=0  THEN  OUTPUT  @Fi lei : Bop . Told* .Emf <*) ,Vr . Ir 

4930  IF  Im=1  OR  Ok  -  I  THEN  OUTPUT  QPlot : Qdp , Thetab 

4940  IF  lm=0  THEN 

4950  BEEP 

4960  INPUT  "WILL  THERE  BE  ANOTHER  RUN  ( 1 =Y . 0=N>?" .Go_on 

4970  Nrun=J 

4980  IF  Go_on<>1  THEM  5030 

4990  IF  Go_on=1  THEN  Repeat 

5000  ELSE 

5010  IF  J<Mrun+1  THEM  Repeat 

5020  END  IF 

5030  IF  Im-0  THEN 

5040  BEEP 

5050  PRINT 

5060  PRINT  USING  "t  OX .  ""MOTE :  "",11,""    data  runs  were  stored  in  file  "".!0A":.J- 

1 .D2_f lleS 

5070  ASSIGN  »File!  TO  * 

5080  OUTPUT  @File2:Nrun-1 

5090  ASSIGN  @File1  TO  D1  file* 

5100  ENTER  <?F  i  lei  ;  DateS  ,Ldtc  1  ,Ldtc2  .  1 1 1 

5M0  OUTPUT  9File2:DateS,Ldtc1  ,Ldtc2.Itt 

5120  FOR  1=1  TO  Nrun-1 

5130  ENTER  ®Fi le I :Bop . Told* ,Emf < » >  ,  Vr  .  Ir 

5140  OUTPUT  eFile2:Bop.Told*.Emf(»).Vr,Ir 

5150  NEXT  I 

5160  ASSIGN  @Fi  le1  TO  * 

5(70  PURGE  "DUMMY" 

5180  END  IF 
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j  i  JU  DLcr 

}2UU  PR1N1 

5210  PRINT  USING  " I  OX  .""NOTE :  "",12,""    X-Y  pairs  were  stored  in  plot  data  file 

"".10ft": J-1 .P_f lleS 

5220  ASSIGN  3»File2  TO  * 

5230  ASSIGN  ©Plot  TO  * 

5240  BEEP 

5250  INPUT    "LIKE    TO   PLOT    Dftfft    ( 1 = Y , 0=N)?M ,0k 

52G0  IF    Ob  =  1    THEN 

5270  CALL  Plot 

5280  END  IF 

5290  SUBEND 

5300! 

5310!  CURVE  FITS  OF  PROPERTY  FUNCTIONS 

5320  DEF  FNKca(T) 

5330!  0FHC   COPPER   250    TO    300   K 

5340  Tk-T+273.15  !C    TO   K 

5350  K=434- . I 12*Tk 

5360  RETURN   K 

5370  FNEND 

5380  DEF  FNMu(T> 

5330!  170  TO  360  K  CURVE  FIT  OF  VISCOUSITY 

5400  Tk  =  T+273 . 15     !C  TO  K 

5410  Mu=EXP<-4.4636+( 1 0  I  1 . 47 /Tk ) )* 1 . 0E-3 

5420  RETURN  Mu 

5430  FNEND 

5440  DEF  FNCp(T) 

5450!  130  TO  400  K  CURVE  FIT  OF  CP 

5460  Tk=T+273.15      !C  TO  K 

5470  Cp». 40 188+1 . 6500  7E-3-Tk + 1 .51494E-6*Tk  2-6 . 67853E- ! U-Tk  3 

5480  RETURN  CP 

5490  FNEND 

5500  DEF  FNRho(T) 

5510  Tk=T+273. 15      !C  TO  K 

5520  X-1-(1 .8*Tk/753.95>   !K  TO  R 

5530  Ro= 36. 32+61 . I46414«X  < 1 /3 >  +  1 6 . 4 1 80 1 5*X+ I  7 . 476838*X  .5+1.1 19828*X'2 

5540  RETURN  Ro 

5550  FNEND 

5560  DEF  FNPr(T) 

5570  Pr=FNCp<T)»FNMu<T>/FNK<T> 

5580  RETURN  Pr 

5590  FNEND 

5600  DEF  FNK(T) 

5510!  T<360  K  WITH  T  INC 

5620  K= . 071  - . 000261 *T 

5630  RETURN  K 

5640  FNEND 

5650  DEF  FNTanh(X) 

5660  P=EXP(X) 

5670  Q=EXP(-X) 

5680  Tanh=(P-Q)/<P+Q> 

5690  RETURN  Tanh 

5700  FNEND 

5710  DEF    FMTvsv(V) 

5720  COM    /Cc/   C(7> 

5730  T=C<0) 

5740  FOR    1=1    TO    7 

5750  T=T+C<I)*V    I 

5  760  NEXT    I 

5  770  T  =  T-6.7422934E-2  +  W9.027  7043E-3-T*(-9.3259917E-5)> 
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5780 

5  790 

5800 

5810 

5820 

5830 

5840 

5850 

5860 

5870 

5880 

5890 

5900 

5910! 

5920 

5930 

5940 

5950 

5960 

5970 

5980 

5990 

6000 

6010 

6020 

6030 

6040 

6050 

6060 

6070 

6080 

6090 

6100 

6110 

6120 

6130 

6140 

6150 

6160 

6170 

6130 

6190 

6200 

6210 

6220 

6230 

6240 

6250 

6260 

6270 

6280 

6290 

6300 

6310 

6320 

6330 

6340 

6350 

6360 

6370 

6380 


RETURN  T 

FNEND 

DEF  FNBeta(T) 

Rop=FNRho<  T+.  I ) 

Rom=FNRho<  T-  .  1 > 

Beta=-2/(Rop+Rom)»(Rop 

RETURN  Beta 

FNEND 

DEF  FNTfilm(Tu,Tld) 

Tfllfn=<Tw+Tld>/2 

RETURN  Tfilm 

FNEND 

DEF  FMPsat(Tc) 

0  TO  80  deg  F  CURVE 

Tf-1 ,8-Tc+32 

Pa=5.9'45525  +  Tf »(  .  15 

Pg=Pa-14.7 

IF  Pq>0  THEN 

Psat=Pg 

ELSE 

Psat=Pq*29.92/14.7 

END  IF 

RETURN  Psat 

FNEND 

SUB  Plot 

DIM  C(9) 

INTEGER  Ii 

PRINTER  IS  1 

Idv-1 

BEEP 

INPUT 

BEEP 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

INPUT 

IF  Opo 

BEEP 

INPUT  "SELECT 

END  IF 

PRINTER    IS    705 

IF    IdvOI    THEN 

BEEP 


Rom)/. 2 


FIT    OF    Psat 
352082+ Tf*< 1 .4840963E-3+rf*9.6150671E-6>> 


!     +=PSIG.-=in    Hg 


LIKE  DEFAULT  VALUES  FOR  PLOT  < I = Y , 0=N>?" , Idv 


USING 

USING 

USING 

USING 

USING  " 

Opo 

■3    THEN 


"4X 
"4X 
"4X 
"4X 
"4X 


"Select  Uption: 

"0 

"1 

"2 

"3 


q  versus  delta-T' 
h  versus  deita-T" 

h  versus  q 

h-ratio  versus 


delta-T 


TUBE  DIAMETER  (0-. 75. 1-1.0  IN)". ltd 


INPUT 

BEEP 

INPUT 

BEEP 

INPUT 

BEEP 

INPUT 

ELSE 

Cy  =  3 

IF  Opc=0 

Cx  =  2 

Xm  i  n= 1 

Yrnin=1000 

END  IF 

IF  0Po=1 

Cx  =  2 


ENTER 
ENTER 
ENTER 
ENTER 

THEM 


THEN 


NUMBER 
NUMBER 
MIN 

MIN 


OF 
OF 


CYCLES 
CYCLES 


FOR 
FOR 


X-AXIS" 

Y-AXIS" 


.C) 


Cy 


X-VALUE    (MULTIPLE    OF    10)",Xmin 
Y-VALUE    (MULTIPLE    OF    lOV.Ymin 
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idoJU  Ami  n=  I 

6400  Ynun=100 

6410  END    IF 

6420  IF  0po=2  THEN 

6430  Cx  =  3 

6440  Xmin=1000 

6450  Ynun-100 

6460  END    IF 

6470  IF  Opo=3  THEN 

6480  Cx=2 

6490  Cy=3 

6500  Xmin=1 

6510  Ymin=l 

6520  END  IF 

6530  END  IF 

6540  BEEP 

6550  PRINT  "IN:SP1:IP  2300,1800.8300.6800 

6560  PRINT  "SC  0  . 1 00  ,  0  . 1 00 : TL  2.0:" 

6570  Sfx=IOO/Cx 

6580  Sfy=100/Cy 

6590  PRINT  "PU  0,0  PD" 

6600  Nn=9 

6610  FOR  1=1  TO  Cx+1 

6620  Xat=Xmin-10  (1-1 > 

6630  IF  I-Cx+1  THEN  Nn-1 

6640  FOR  J=1  TO  Nn 

6650  IF  J-1  THEN  PRINT  "TL  2  0" 

6660  IF  J  =  2  THEN  PRINT  "TL  1  0" 

6670  Xa=Xat*J 

6680  X-LGT<Xa/Xmin>»Sfx 

6690  PRINT  "PA":X.",0:  XT;" 

6700  NEXT  J  • 

6710  NEXT  I 

6720  PRINT  "PA  I00,0;PU:" 

6730  PRINT  "PU  PA  0,0  PD" 

6740  Nn  =  9 

6750  FOR  1=1  TO  Cy+1 

6760  Yat=Ymin-10  (1-1 > 

6770  IF  I-Cy+1  THEN  Nn= I 

6780  FOR  J-1  TO  Nn 

6790  IF  J-1  THEN  PRINT  "TL  2  0" 

6800  IF  J  =  2  THEN  PRINT  "TL  1  0" 

6810  Ya=Yat*J 

6820  Y=LGT< Ya/Ym in) *Sf y 

6830  PRINT  "PA  0,":Y."YT" 

6840  NEXT  J 

6850  NEXT  I 

6860  PRINT  "PA  0,  100  TL  0  2" 

6870  Nn=9 

6880  FOR  1=1  TO  Cx+1 

6890  Xa1-  =  Xmin»10  (1-1) 

6^00  IF  I=Cx+!  THEN  Nn-1 

6910  FOR  J-1  TO  Nn 

6920  IF  J-1  THEN  PRINT  "TL  0  2" 

6930  IF  J>1  THEN  PRINT  "TL  0  1" 

6940  Xa=Xat«J 

6950  X=LGT(Xa/Xmin)*Sfx 

6960  PRINT  "PA";X.", 100:  XT" 

6970  NEXT  J 

6930  NEXT  I 
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6990 

7000 

7010 

7020 

7030 

7040 

7050 

7060 

7070 

7080 

7090 

7100 

7!  10 

7120 

7130 

7140 

7150 

7160 

7170 

7180 

7190 

7200 

7210 

7220 

7230 

7240 

7250 

7260 

7270 

7280 

7290 

7  3J0O 

7310 

7320 

7330 

7340 

7350 

7360 

7370 

7380 

7390 

7400 

7410 

7420 

7430 

7440 

7  450 

7460 

7470 

7480 

7490 

7500 

7510 

7520 

7530 

7540 

7550 

7560 

7570 

7580 

7590 


PRINT  "PA  100,100  PU  PA  100.0  PD" 

Nn  =  9 

FOR  1-1  TO  Cy+1 

Yat  =  Yrmn*10  (1-1) 

IF  I-Cy+1  THEN  Nn= I 

FOR  J-1  TO  Nn 

IF  J-1  THEN  PRINT  "TL  0  2" 

IF  J>1  THEN  PRINT  "TL  0  1" 

Ya=Yat*J 

Y-LGT<Ya/Ymin)-Sfy 


PRINT 
NEXT  J 
NEXT  I 
PRINT 
PRINT 


PD  PA  100 


Y  "YT' 


"PA  100,100 
"PA  0,-2  SR 
Ii-LGT(Xmin) 
FOR  1=1  TO  Cx+1 
Xa=Xmin»!0  (1-1  ) 
X=LGT(Xa/Xmin)»Sfx 
PRINT  "PA";X,",0:" 
IF  li>=0  THEN  PRINT 
IF  li<0  THEN  PRINT 
Ii-H  +  1 
NEXT  I 

PRINT  "PU  PA  0,0" 
Ii=LGT(Ymin> 
Y10-10 

FOR  1=1  TO  Cy+1 
Ya=Ynun«!0  (1-1) 
Y=LGT(Ya/Ymin)*Sfy 


PU' 


1  .5 


"CP 

•cp  ■ 


-2;LB10:PR  -2 
-2:LB10:PR  0.2: 


2;LB";Ii 
LB"  :  1 1  :  "' 


'PA 
CP 


PRINT 
PRINT 
Il-Ii+1 

NEXT  I 
IF  IdvOI 
BEEP 
INPUT 
BEEP 
INPUT 
ELSE 

IF  Opo=0 
XlabelS= 
YlabelS= 
END  IF 
IF  0po=1 
XlabelS= 
YlabelS= 
END  IF 
IF  0Po=2 
XlabelS= 
YlabeU- 
END  IF 
IF  Opo-3 
XlabelS= 
YlabelS= 
END  IF 
END  IF 
PRINT  "SR 
PRINT  "PA 
PRINT  "CP 
Repeat : ! 


0 


■4.-.25;LB10;PR  -2.2;LB":Ii 


THEN 


"ENTER  X-LABEL".Xlabel$ 
"ENTER  Y-LABEL".YlabelS 


THEN 

Tw-Tsat  (KV 
q  (W/m  2)" 

THEN 

Tu-Tsat  <K)' 
h  (W/m  2.K)' 


THEN 
q  (W/m 
h  (W/m 


2)" 
2.K)' 


THEN 

Tu-Tsat  <K>" 
h(enh) /h( smooth)' 


!  .5,2:PU  PA 
-14.50  CP  0 
0.0  DI" 


50. -16  CP": -LEN(X labelS) /2:"0: LB"; X  labels 
":-LEN(YlabelS)/2»5/6;"DI  0.1 : LB" ; Y label? 
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/bUU 

7610 
7620 
7630 
7640 
7650 
7660 
76  70 
7G80 
7690 
7700 
7710 
7720 
7730 
7740 
7750 
7760 
7770 
7780 
7790 
7300 
7810 
7820 
7830 
7840 
7850 
7860 
7370 
7880 
7890 
7900 
7910 
7920 
7930 
7940 
7950 
7960 
7970 
7980 
7990 
3000 
8010 
3020 
3030 
3040 
3050 
3060 
80  70 
3030 
3090 
8100 
31  10 
8120 
3130 
8140 
3!  50 
3160 
8170 
3130 
8190 


0  PLOT  DATA  FROM  Pi  FILE  (1=Y.0=N)?' 


Ok 


THE    NAME    OF 
TO    D    fileS 


THE    DATA    FILE".D    file? 


tittf 

INPUT    "WANT    1 

IF    Dlc-1    THEN 

BEEP 

INPUT    "ENTER 

ASSIGN  SFile 

BEEP 

BEEP 

INPUT  "ENTER  THE  BEGINNING  RUN  NUMBER", Md 

BEEP 

INPUT  "ENTER  THE  NUMBER  OF  X-Y  PAIRS  STORED" . Npai rs 

BEEP 

PRINTER  IS  1 


PRINT  USING  "4X.  ""Select 

PRINT  USING  "4X.""1 

PRINT  USING  "4X.""3 

PRINT  USING  "4X.""5 

PRINT  USING  "4X.""6 

PRINT  USING  "4X.""7 

INPUT  Sym 

PRINTER  IS  705 

PRINT  "PU  DI" 

IF  Sym=1  THEN  PRINT 

IF  Syrn  =  2  THEN  PRINT 

IF  Sym=3  THEN  PRINT 

IF  Md> 1  THEN 

FOR  1=1  TO  (Md-1) 

ENTER  «>File:Ya.Xa 

NEXT  I 

END  IF 

FOR  1=1  TO  Npairs 

ENTER  °>File:Ya.Xa 

IF  0po=1  THEN  Ya=Ya/Xa 

IF  0po=2  THEN 

0  =  Ya 

Ya=Ya/Xa 

Xa=Q 

END  IF 

IF  0Po=3 

X 

Y 

IF  Sym>3  THEN  PRINT  "SM" 

IF  Sym<4  THEN  PRINT  "SR 

PRINT  "PA".X.Y."" 

IF  Sym>3  THEN  PRINT 

THEN 

THEN 

THEN 

THEN 


svmbo  1  : 

Star      2   Plus  sign 

Circle    4   Square 

Rombus 

R  iqht-side-up    tnangl 
Up-side-doun    triangle 


"SM»" 
"SM+" 
"SMo" 


e 

IIKII 


THEN    Ya=Ya/FNHsmooth(Xa.Itd) 
LGT(Xa/Xmin)»Sfx 
LGT(Ya/Ymin)»Sfy 


1  .4 


4" 


IF 

IF 
IF 
IF 


PRINT 
PRINT 
PRINT 
PRINT 


"SR  1 
"UC2. 
"UC3, 
"UC0  . 
"UC0, 


99 
99 


.6" 

,o.-t 

.  -  3  ,  - 

99.3 
3.99.- 


,-4 

6,-: 

-8.- 

3,S 


0.0 
.6, 
6.0 
6.0 


Sym  =  4 

Sym  =  5 

Sym  =  6 

Sym=  7 
NEXT  1 
BEEP 

ASSIGN  ©File  TO  * 
GOTO  7600 
END  IF 

PRINT  "PU  SM" 
BEEP 

INPUT  "WANT  TO  PLOT 
IF  Go_on=1  THEN 
BEEP 
INPUT  "ENTER  LOWER  AND  UPPER  X-LIMITS" . XI 


8.4.0 

i. 6 . 3  .  - 

3.8:" 

-3,-8 


A    POLYNOMIAL    < I = Y , 0=N) ?" , Go. 

Xlu 
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8200  FOR  Xx  =  0  TO  Cx  STEP  Cx/200 

8210  Xa=Xmin*!0  Xx 

8220  IF  Xa<Xll  OR  Xa>Xlu  THEN  8290 

8230  Ya=FNPoly<Xa) 

8240  Y=LGT<Ya/Ymin)»Sfy 

8250  X=LGT(Xa/Xmin)*Sfx 

8260  IF  Y<0  THEM  Y=0 

3270  IF  Y > 1 0 0  THEN  GOTO  8290 

8280  PRINT  "PA'\X,Y,"PD" 

8290  NEXT  Xx 

8300  END  IF 

8310  PRINT  "PU  PA  0.0  SPO" 

8320  SUBEND 

3330  DEF  FNHsmooth(X.Itd) 

3340  Hs=FNPoly(X)/X 

3350  IF  ltd- 1  THEN  Hs=Hs*. 83347 

8360  RETURN  Hs 

3370  FNEND 

3380  DEF  FNPoly(X) 

3390  Poly=-4.4123713E+2-X»(6.8123917E+2-X»3.7416863E+2> 

8400  RETURN  Poly 

3410  FNEND 

9420  DEF  FNPvst(Tsteam) 

8430  DIM  K(8> 

3440  DATA  -7.691234564,-26.08023696,-168. 1706546.64.23285504,-1 18.9646225 

8450  DATA  4 . 1 671 1 732 .20 . 9750676  .  1  . E9  .  6 

3460  READ  K(-) 

3470  T=<Tsteam+273.15)/647.3 

8480  Sum=0 

8490  FOR  N-0  TO  4 

3500  Su.m  =  Sum+K<N)*< 1  - T  >  C N  + 1  ) 

8510  NEXT  N 

3520  Br  =  SLLm/(T»(  1+K(5)*(  1-T)+K<6)*(  l-T)    2  >  )  -  (  1  -  T  )  /<  K  <  7  )  *(  1  -T  )    2+K<8)> 

8530  Pr=EXP(Br) 

3540  P=22120000«Pr 

8550  RETURN   P 

3560  FNEND 
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APPENDIX  F 
AN  EXAMPLE  OF  REPRESENTATIVE  DATA  RUN 


Month,  date  and  time  : 12 : 03: 16 : 50 :24 

NOTE:  Program  name  :  DRP 
Disk  number   *  02 
Old  file  name:  NH05 

This  data  set  tak»n  on  :  11:02:10:36:28 
Thermocouples  were  defective  at  locations:   0       0 
Tube  Type  l s  :  2 

Oat  a  Set  Number  -    !   Bulk  Oil  '/.    -   0.0      I  I  :02 :  1  t  :  0 1  :46 

TC   No:  1  2  3  4  5  6  /  B 

Temp    :       70.  OS      S5.ll  7      (51.04      7U.?3      78.68      SO..:.?      78.21      59.24 

Twa  find        Thetab  Tvapr      Pea!'         Isump 

66.34      10.03      5.560E+0  1         10.01         3.86      -13.0 

Htube  Qdp 

I.545E+03      8.588E+04 

Oata   Set    Number    -        2      Bulk    Oil    '/.    =      0.0  11:02:11:30:52 

TC    No:  1  2  3  4  5  6  7  F: 

Temp    :      52.64      42.70      45.49      49.95      58.26      44.85      57.78      47.04 
Twa  TUad         Thetab  Tvapr      Psat         Tsump 

49. 84       1 0.0 7      3.939E+01         10.30        3.88      -12.0 
Htube  Qdp 

.1 . I52E+03      4.538E+04 

Data   Set    Number    ■        3      Bulk    Oil    7.    -      0.0  11:02:11:38:07 

TC    No:  12  3  4  5  6  7  8 

Tenp    :      34.63      30.25      31.40      33.49      37.59      3  1,38      38.00      32.96 

Twa     Tliqd    Thetab      Tvapr   Psat    Tsump 

3  3.72   '0.00   2.357E  +  01    10.69    3.8:3   -11.9 

Htube      Qdp 

7.420E+02   1 .749E+04 

Data  Set  Number  -    4   Bulk  Oil  7.  -   0.0      11:02:11:49:56 

TC  No:     1       2      3      4  f:.      7      8 

Temp  :   27.30   24.92   25.40   26.60   29.04   25.25   29.30   26.40 

Iwa     T 1 i qd    Thetab      Tvapr   Psat    Tsump 

26.78   10.07   1.663E+01    11.17   3.88   -11.9 

Htube      Qdp 

5.2HE  +  02     8.664E+03 

Data   Set    Number    -         5      Bulk    Oil    X    -      0.0    .  1  '  : 02  :  1  I  :  57  :  02 

TC   No:  1  2  3  4  5  fi  7  8 

Tenp    :      22.10      21.03      21.06      21.41      23.11      21.55      23.11      21.65 

Twa     T 1  iqd    Thetab      Tvapr   Psat    Tsump 

21.99   9.95   1.189E+01    I '.4 7   3.80   -12.1 

Htube       Qdp 

3.094E+02   3.630E+n3 

Data  Set  Number  -    6   Bulk  Oil  V.    -   0.0      11:02:12:05:48 
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TC  No:     1      2      3      4      5      6      7      ft 
Temp  :   18.14   17.58   17.82   17.57   19.16   '8.03   18.43   16.84 

Twa     Tliqd    The  tab      Tvapr   Psat    Tsump 

17.95   10.19   7.737E+00    12.00    3.36   -12.2 

Htu.be       Odp 

I.728E+02      I.337E+03 


Data   Set    'lumber 
1 


7   Bulk  Oil  /. 
3 


0.0 


02:12: 14: 18 


TC  No:     1       2      3      4      5      6      7 

Temp  :  1 5 . 58   1 5 .  08   15.14   14.37   15.86   15.04   15.46 

Twa  Tliqd    Thetab      Tvapr   Psat    Tsump 

15.17  10.12  5.042E+00    12.34   3.92   -12.2 

Htube      Odp 

I.221E+02  6. I56E+02 

NOTE:  07  X-f  pairs  mere  stored  in  plot  data  file  P05 
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too 

1  10 

120 

130 

140 

ISO 

160 

170 

180 

190 

200 

21  0 

220 

230 

2^0 

241 

242 

244 

245 

250 

260 

270 

280 

290 

300 

310 

320 

3^0 

340 

350 

360 

361 

362 

363 

364 

365 

3  70 

3R0 

390 

4  00 

410 

411 

4I2 

413 

420 

430 

440 

450 

4  60 

470 

480 

490 

500 

510 

520 

530 

5  40 

541 

542 

544 

APPENDIX    G 
LISTING    OF    CALIBRATION    COMPUTER    PROGRAM    (TCAL) 

•    FILE    NAME:    TCAL 

!    REVISED:         October    4.     198^ 
I 

COM    /Cc/    C<7) 
DIM    EmF(4) , T(4) ,D(4) 

DA  I  A  0.  10  086091  .  25  72  7  .  94  369  ,-  76  7345  .  8295  .  78025585  .  8  1 
DATA  -924 7486589. 6. 97688E1 1  , -2 . 661 92E 1 3 . 3 . 940  78E 1 4 
READ  C<*> 
CLEAR  709 
BEEP 

INPUT  "ENTER  MflNTH.  DATE  AND  TIME  < MM : DD : HH : MM : SS ) " , BS 
J  =  0 

OUTPUT  709;"TD":B$ 
OUTPUT  709:"TD" 
ENTER  709;A$ 
BEEP 

INPUT  "WANT  A  HARD  COPY  <  1  -  Y  ,0-N)?" , Ihp 
IF  Ihp-1  THEN  PRINTER  IS  701 
IF  Ihp-0  THEN  PRINTER  IS  1 

PRINT  USING  "10X, ""Month,  date  and  time   =  "".14A";A$ 
BEEP 

INPUT  "ENTER  INPUT  MODE  (1=3054A.  2=FILE>" . Im 
IF  Im-1  THEN 
BEEP 

INPUT  "GIVE  A  NAME  FOR  DATA  F ILE" , D_f i leS 
CREATE  BOAT  D_fileS,20 
ELSE 
BEEP 

INPUT  "GIVE  NAME  OF  EXISTING  F ILE" , D_f i 1 e$ 
BEEP 

INFUT  "ENTER  NUMBER  OF  DATA  RUNS  STORED", Nrun 
BEEP 

INPUT  "GIVE  A  NAME  FOR  OUTPUT  FILE"  .Of i  leS 
CREATE  BDAT  Of  UpS,  10 
ASSIGN  ?Fileo  TO  Of lle$ 
K  =  1 

END  IF 

ASSIGN  QV i le  TO  D_f  i le« 
IF  Im-1  THEN 
BEEP 
!  INPUT  "ENTER  BATH  TEMPERATURE" , T_bath 
OUTPUT  713;"T2R2E" 
WAIT  2 

ENTER  713;T_bath 
OUTPUT  709:"AR  AF29  AL 33  VR 1" 
FOR  1=0  IN  4 
OUTPUT  70q:"AS  SA" 
ENTER  709;Emf(D 
NEXT  I 
ELSE 

ENTER  •'File:T_bath,Emf<-> 
END  IF 
J-J+1 

FOR  1=0  TU  4 

T<  I)=FNT>/?v<ABS<Emf  (  I)  )) 
D(I)=T_bath-T<I) 
NEXT    I 
IF    K=1     THEN 
FOR    L=0    TO    4 
OUTPUT    flFileo:T(L) ,D<L) 
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545 
548 
549 

550 

560 
570 
530 
530 
600 
610 
620 
630 
640 
650 
660 
670 
630 
630 
700 
710 
720 
730 
740 
750 
760 

7  70 
780 
730 
_f  i 
300 
810 

f  i  le 

8  1  1 
820 
330 
840 
350 
860 
3  70 
880 
R^O 
900 
310 
920 


NEXT  L 

GOTO  700 

END  IF 

PRINT 

PRINT  USING  "10X.""Bath  temperature       -  ""  .3D. 3D 

PRINT  USING  "1  OX . ""Thermocouple  readinqs  <Deg  C>: 

PRINT  USING  "inx.5<3D.DD,3X> ,25X";T<«> 

PRINT  USING  "10X, ""Discrepancies  <Deg  C): 

PRINT  USING  "1  IX.5(r1Z.DD.4X>  .24X";D(»> 

IF  Im=2  THEN  700 

BEEP 

INPUT  "UK  TO  ACCEPT  THIS  SET  ( 1 -Y,0-N)?".0ks 

IF    Dlfs-O    THEN 

J-J-1 

GOTO    390 

ELSE 

OUTPUT  $Fi le; T_bat h . Emf < * ) 

END  IF 

BEEP 

IF  Im=1  THEN 

INPUT  '■WILE  THERE  BE  ANOTHER  SET  <  1  -Y  ,  0-N)?M  ,Go_on 

IF  Go_on=1  THEM  390 

ELSE 

IF  J<Nrun  THEN  390 

END  IF 

PRINT 

IF  Im=1  THEM 

PRINT  USING  ' 


(Deq  C) ; T_bath 


le$ 


ELSE 

PRINT  USING 


?S 


OX 

OX 

ox 


MOTE 
NOTE 


DD. 


data  sets  are  stored  in  file  "".I4A";J,D 


PRINT  USING 

END  IF" 

ASSIGN  fr'i  le'TQ  « 

END 

DEF  FNTvsv(Fmf ) 

COM  /Cc/  C(7) 

Sum'C(O) 

FOR  1-1  TO  7 

Sum  =  Sum«-C(I)  "Emf  I 

NEXT  I 

RETURN  Sum 

FNEND 


Above  analysis  uas  performed  from  file  "",!4A";D_ 
Output  data  are  stored  in  file  "" . ! 4A" ; Of i leS 
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APPENDIX  H 
SAMPLE  CALCULATION 

Data  run  number  5  (saturation  temperature  was  10  °C  and 
heat  flux  was  about  20  kW/m2 )  was  chosen  for  the  sample 
calculation . 

A.   TEST- SECTION  DIMENSIONS 


Di 

L 
L 


u 


0.01588  (m) 

0.01270  (m) 

0.01143  (m) 

0.01588  (m) 

0.20320  (m) 

0.07620  (m) 


B.   MEASURED  PARAMETERS 


vs 

= 

3.58    ( 

^volts ) 

Is 

= 

2.27    ( 

Volts ) 

Ti 

- 

34.69    ( 

;°c) 

T2 

= 

30.25    ( 

;°c) 

T3 

- 

31.40    ( 

;°c) 

T4 

- 

33.49    ( 

;°c) 

T5 

- 

37.59    ( 

:°c) 

T6 

- 

31.38    1 

;°c) 

T7 

- 

38.00    1 

;°c) 

T8 

- 

32.96 

;°c) 

T  sat 

= 

10.00    1 

;°c) 

kc 

= 

344.00 

;W/m.K) 
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C.  OUTER  WALL  TEMPERATURE  OF  THE  BOILING  TUBE 

p   =  *    ^   =  3.1416  x  0.01588  =  0.0499  (m) 

Ac  =  *     (D20  -    D2±)/4    =    3.1416    x    (0.015882    -    0.01272)/4 

Ac  =       714    x    10  -7     (m2) 

V      =      25VS  =    25    x    3.58    =    89.50    (volts) 

I      =      Ig  =    2.27    (amp) 

Q      =      VI    =    89.50    x   2.27    =    203.17    (W) 

8 

^vg=       (    £  TLJ    /    8    =    33.72    (°C) 
n=i 

f      =      T    -    Q    (    ln(D  J    DJ/2   ^L   k    ) 
wo  avg  XHV  K     2'         V  c  ' 

\0-      33.72    -    203. 17(ln(15. 88/11. 43)/2x3. 1416x203.2x344) 
T      =       33.57    (°C) 

wo  v         ' 

Tb     "       Two 

9      =      Tb-    T      =    3.3.57    -    10.00    =    23.57    (°C) 
"        sat 

D.  PROPERTIES    OF    R-114   AT    FILM   TEMPERATURE    {REF.    24},     {REF. 
25} 

T-  =       (f  +    T     J/2    =    21.78    (°C)       =      294.78    (K) 

f  w  o       s  a  tr 

log  y (10  "3Ns/m2)    =    -4.4636    +    1011.47/Tf 

y       =    exp(-4.4636    +    (1011.47/294.78))    x    10  "3 

y       =    356    x    10  "6     (Ns/m2) 

P       =       36.32    +    61.14   j1/3    +    16.42    j    +    17.48    j1'2    + 

1.12    j2       (lb/ft3) 
J       --       1    -    Tf  /T 

T   =   Critical  Temperature  =  753.95  (R) 
T   =   294.78  (K)  =   530.6  (R) 
j   =   1  -  530.60/753.95  =  0.296 
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91.55  (lb/ft3) 


1466.49  (kg/m3) 


v   =   U  /  P  =   2.428  x  10  "7  (m2/s) 

k  (W/m.K)   =  0.0710  -  0.000261  Tf    (Tf  in  °C) 
k   =   0.0653  (W/m.K) 


C    (kJ/kg.K) 


=    0.4    +    1.65    x    10"3Tf+    1.51    x    10  -eT. 
6.68    x    10  -1  °Tf3       (T    in  K) 


(^  =       1002    (J/kg.K) 

a  =      k/  P  Cp  =    4.444   x    10  "8     (m2/s) 

6  =  -  (A  p  /  A  T )  /   P 

6  =  -(1465.95    -    1466.49)    /    ((294.98    -    294.78)    x    1466) 

B  =      0.00184    (1/K) 


Pr 


v /  a    =      5.464 


E.   HEAT -FLUX  CALCULATION 

Average  natural- convection  heat- transfer  coefficient 
at  non-boiling  ends: 


h      = 


D 


<  0.60  +0.387 


gSD3eRTanh;(Jl£-)1/2Lut11/6 

^     O   D  (  kQ^Q 

vaLu(r-A-} 

C  C 


[1  +  (0.559/Pr)9/16]8/27 


h   =   282  (W/m2.K) 

Heat- transfer  rate  through  non-boiling  ends 
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(hpkA    )"2    9,Tanh(mL    ) 


m 


(h   p    /    kcAc) 


1     /  2 


23.94    (1/m) 


0      =      13.15    (W) 
0      _=      2    x   0^=    2    x    13.15      =    26.30    (W) 


Heat  flux  through  active  boiling  surface 


Q 
q 

A* 
q 

h 


0    -    Q^      =    203.17    -    26.30    =    176.87    (W) 
Q   /   Ab 
D   L    =    10.  137    x    10  "3     (m2) 


176.87  /  0.010137 


17447.96  (W/m2) 


q  /  ev 


740.26  (W/m2.K) 


The  following  are  the  results  obtained  from  the  computer 
by  running  the  data  reduction  program  (See  Appendix  F). 


9, 


17490  (W/m2) 
23.57  (°C) 
742  (W/m2.K) 
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APPENDIX  I 
UNCERTAINTY  ANALYSIS 

The  same  data  set  (run  number  5)  that  was  used  for  the 
sample  calculation  was  chosen  for  the  uncertainty  analysis; 
therefore,  the  dimensions  of  the  test  section,  and  the  meas- 
ured and  calculated  parameters  found  in  the  sample  calcula- 
tion were  used  in  this  analysis.  All  uncertainties  are 
presented  as  a  percentage  of  the  calculated  parameter. 

A.   UNCERTAINTY  IN  SOURCE  HEAT-TRANSFER  RATE 


QH  =   VI  (W) 

I   =   I   =  2.27  (amp)  <$I   =   ±0.025  (amp) 

V   =   25V  ;  V   =  3.58  (volts)    6VS  =   ±0.05  (volts) 
s    s 

where : 

5   =   uncertainty  in  measurement  and  calculation 

6QH/QH  =   ((5VS/VS)2  +  (Sl/l)2)i  " 

5QH/QH  =   ((0.05/3.58)2  +  (0 . 025/2 . 27  ) 2 )  l  '  2 
-      1.80  percent 

B.   UNCERTAINTY  IN  SURFACE  AREA 


Ab    =       uDQL 

DQ    =         15.88     (mm)  6DQ   =       0.1    (mm) 

L      =       203.20    (mm)  6L      =       0.1    (mm) 

6Ab/Ab    =       ((5Do/Do)2    +     (5L/L)2)1'2 

5Ab/Ab    -       ((0.1/15.88)2    +    (0. 1/203. 2)2)1  '2 
=       0.63    percent 


111 


UNCERTAINTY  IN  WALL  SUPERHEAT 


AT   =   T^0  -   T 


sat 


T      =       10.00    (°C)  £T      =      0.5    (°C) 

^o  =      W    Qh  (InCPg  /Dx  )    /    2  tt  L   kc  ) 


t\) 


T       =          )    T    J   /    8 
avg      \  L i 

n=  1 

where : 

T^  =   thermocouple  readings  (See  Appendix  G) 
Tav#   33.72  («C) 

S.D.  =   ((^V-  Tavg)2)/7  )1/2   =   2.687   (°C) 

n=l 
where : 

S.D.  =   standard  deviation 

Since  logarithmic  term  in  equation  of   "T^  "  is  too  small 

when   compared  to   standard   deviation,    this  term   can   be 

neglected  for  uncertainty  analysis  (i.e.,  TL  =   T    ). 

T   =   33.72  (°C)  <5f   =   2  x  S.D.  =  5.37  (°C) 

wo  wO 

aT      =       f       -       T         =       33.72    -    10.00    =    23.72    (°C) 

u  wo  sat  v       ' 

6AT/AT      =       ((6WAT)2    +     (  -5  T^/AT  ) 2  )  l  '  2 

5AT/AT      =       ((5. 37/23. 72)2    +    ( -0 . 5/23  .  72  ) 2  )  x  '  2 
6AT/AT      =       22.7    percent 

D.       UNCERTAINTY    IN   HEAT    FLUX 


Q      =      203.17    (W)  6Q      =      3.65    (W) 

^  H 

Assuming    the    same    proportion    in    the    uncertanity    for    Q: 
Qp     =       13.15    (W)  6Q      =      0.23    (W) 

Qj     -      2QF    =       176.87    (W) 
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6q/q   =  (  (5QR/(QH-  2Q_f))2  ♦  (26Qp/(QH-  2Qp  ))2  + 
(5i^/^)2  )»  '2 

6q/q   =  (  (3. 65/176. 87)2  +  (0 . 46/ 176 . 87 ) 2  + 
(-0.0063)2  )'  '2 
=       2.17  percent 

E.   UNCERTANITY  IN  BOILING  HEAT-TRANSFER  COEFFICIENT 

h   =   q  /  A  T 
fih  /  h   =   (  (6q/q)2  +  (-5AT/AT)2  )»  '2 

5h  /  h   =   (  (0.0217)2  +  (-0.227)2  )1/2 
=   22.8  percent 
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